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Abstract 

In non-relativistic mechanics the center of mass of an isolated system is easily separated out 
from the relative variables. For a N-body system these latter are usually described by a set of 
Jacobi normal coordinates, based on the clustering of the centers of mass of sub-clusters. The 
Jacobi variables are then the starting point for separating orientational variables, connected with 
the angular momentum constants of motion, from shape (or vibrational) variables. Jacobi variables, 
however, cannot be extended to special relativity. We show by group-theoretical methods that two 
new sets of relative variables can be defined in terms of a clustering of the angular momenta of 
sub- clusters and directly related to the so-called dynamical body frames and canonical spin bases. 
The underlying group-theoretical structure allows a direct extension of such notions from a non- 
relativistic to a special-relativistic context if one exploits the rest-frame instant form of dynamics. 
The various known definitions of relativistic center of mass are recovered. The separation of 
suitable relative variables from the so-called canonical internal center of mass leads to the correct 
kinematical framework for the relativistic theory of the orbits for a N-body system with action-at- 
a-distance interactions. The rest-frame instant form is also shown to be the correct kinematical 
framework for introducing the Dixon multi-poles for closed and open N-body systems, as well as 
for continuous systems, exemplified here by the configurations of the Klein-Gordon field that are 
compatible with the previous notions of center of mass. 

Invited contribution for the book Atomic and Molecular Clusters: New Research (Nova Science). 
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I. INTRODUCTION. 



In nuclear, atomic and molecular physics, as well as in celestial mechanics, an old basic 
problem is to exploit the translational and rotational invariance of a N-body system for 
eliminating as many global variables as possible to get a reduced system described by a 
well defined set of relative degrees of freedom. In molecular dynamics, for instance, this 
reduction is instrumental for the definition of molecular vibrations and rotations. A similar 
old problem exists for isolated continuous deformable bodies: in this case the usual Euler 
kinematics of rigid bodies in the body frame must be generalized to understand phenomena 
like the falling cat and the diver. We suggest to see the review part of Ref. [1] for an essential 
selection of the huge bibliography on such issues. 

Since most of the applications concerning these problems are typically non-relativistic, 
nearly all treatments deal with non-relativistic systems of the kinetic-plus-potential type. 
However, developments in particle physics, astrophysics and general relativity suggest to 
extend the treatment from the absolute Galilei space plus time to Minkowski space-time 
and then to Einstein space-times. Furthermore, the case of the general relativistic N-body 
problem, in particular, forces us to simulate the issue of the divergences in the self-energies 
with a multipolar expansion. We need, therefore, a translation of every result in the language 
of such kind of expansions. 

In this article we show a new method for the treatment of relative variables for a many- 
body system both non-relativistic and relativistic. Our proposal is based upon a conjunction 
of Hamiltonian and group-theoretical methods leading to a systematic generalization, valid 
for generic deformable systems, of the standard concept of body frame for rigid systems. We 
also show that our procedure constitutes the proper kinematical framework for introducing 
multipolar expansions for closed and open systems. 

In Newton mechanics isolated systems of N particles possess 3N degrees of freedom in 
configuration space and 6N in phase space. The Abelian nature of the overall translational 
invariance, with its associated three commuting Noether constants of the motion, makes pos- 
sible the decoupling and, therefore, the elimination of either three configurational variables 
or three pairs of canonical variables, respectively (separation of the center- of-mass motion). 
In this way one is left with either 3N-3 relative coordinates p a or 6N-6 relative phase space 
variables p a , p a , a = 1, .., N — 1 while the center-of-mass angular momentum or spin is 
S = a=i Pa x Pa- Most of the calculations of the non-relativistic theory employ the sets 
of 3N — 3 Jacobi normal relative coordinates s a (see Ref.[l]) that diagonalize the quadratic 
form associated with the relative kinetic energy (the spin becomes S = Yla=i x ^sa, with 
n sa momenta conjugated to the s a 's). Each set of Jacobi normal coordinates s a orders the 
N particles into a hierarchy of clusters, in which each cluster of two or more particles has 
a mass given by an eigenvalue (reduced mass) of the quadratic form; the Jacobi normal 
coordinates join the centers of mass of cluster pairs. 

On the other hand, the non- Abelian nature of the overall rotational invariance entails that 
an analogous intrinsic separation of rotational (or orientational) configurational variables 
from others which could be called shape or vibrational^ impossible. As a matter of fact, this 
is one of the main concerns of molecular physics and of advanced mechanics of deformable 
bodies. In fact, in the theory of deformable bodies one looses any intrinsic notion of body 
frame, which is the fundamental tool for the description of rigid bodies and their associated 
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Euler equations. A priori, given any configuration of a non-relativistic continuous body (in 
particular a N-body system), any barycentric orthogonal frame could be named body frame 
of the system with its associated notion of vibrations. 

This state of affairs suggested [1] to replace the kinematically accessible region of the non- 
singular configurations [2] in the (3N-3)-dimensional relative configuration space by a SO (3) 
principal fiber bundle over a (3N-6)-dimensional base manifold, called shape space. The 
SO (3) fiber on each shape configuration carries the orientational variables (e.g. the usual 
Euler angles) referred to the chosen body frame. Then, a local cross section of the principal 
fiber bundle selects just one orientation of a generic N-body configuration in each fiber 
(SO (3) orbit). This is in fact equivalent to a gauge convention namely, after a preliminary 
choice of the shape variables, to a possible definition of a body frame {reference orientation). 
It turns out that this principal bundle is trivial only for N=3, so that in this case global cross 
sections exist, and in particular the identity cross section may be identified with the space 
frame. Any global cross section is a copy of the 3-body shape space and its coordinatization 
provides a description of the internal vibrational motions associated with the chosen gauge 
convention for the reference orientation. For iV > 4, however, global cross sections do 
not exist [4] and the definition of the reference orientation {body frame) can be given only 
locally. This means that the shape space cannot be identified with a (3N -6) -dimensional 
sub-manifold of the (3N-3)-dimensional relative configuration space. The gauge convention 
about the reference orientation and the consequent individuation of the internal vibrational 
degrees of freedom requires the choice of a connection T on the SO (3) principal bundle (i.e. a 
concept of horizontality) , which leads in turn to the introduction of a SO (3) gauge potential 
on the base manifold. In this way a natural gauge invariant concept exists of purely rotational 
N-body configurations {vertical velocity vector field, i.e. null shape velocities). Clearly, a 
gauge fixing is needed in order to select a particular local cross section and the correlated 
gauge potential on the shape space. Obviously, physical quantities like the rotational or 
vibrational kinetic energies and, in general, any observable feature of the system must be 
gauge invariant. On the other hand, in the orientation-shape bundle approach, both the 
space frame and the body frame components of the angular velocity are gauge quantities 
so that their definition depends upon the gauge convention. See Ref.fl] for a review of the 
gauge fixings used in molecular physics' literature and, in particular, for the virtues of a 
special connection C corresponding to the shape configurations with vanishing center-of- 
mass angular momentum S. The C connection is defined by introducing local cross sections 
orthogonal to the fibers with respect to the Riemannian metric dictated by the kinetic 
energy. 

The orientation-shape approach replaces the usual Euler kinematics of rigid bodies and 
implies in general a coupling between the internal shape variables and some of the orienta- 
tional degrees of freedom. In Ref. [1] it is interestingly shown that the non-triviality of the 
SO (3) principal bundle, when extended to continuous deformable bodies, is the core of the 
physical explanation of the falling cat and the diver. A characteristic role of SO (3) gauge 
potentials in this case is to generate rotations by changing the shape. This approach is de- 
fined in configuration space and leads to momentum-independent shape variables, so that it 
can be easily extended to the quantum level (the N-body Schrodinger equation). In Ref. [1] 
the Hamiltonian formulation of this framework is also given, but no explicit procedure is 
worked out for the construction of a canonical Darboux basis for the orientational and shape 
variables. See Refs.fl, 3] for the existing sets of shape variables for N = 3,4 and for the 
determination of their physical domain. 
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Consider now the description of a N-body system in a special relativistic context. Unlike 
the Newtonian case, manifest Lorentz covariance requires the introduction of 4N degrees of 
freedom in configuration space and 8N in phase space. In phase space there are N mass- 
shell first-class constraints (epj — mf ps in the free case [5]). They determine the energies 
p° ± \/mf + pf and entail that the time variables x° be gauge variables. Such gauge 
variables can be replaced by a center- of-mass time (arbitrariness in the choice of the rate of 
the center-of-mass clock) and N — 1 relative times (arbitrariness in the synchronization of 
the clocks associated to each body). The problem of relative times has been a big obstacle to 
the development of relativistic mechanics (see Refs.[6] for the essential bibliography). The 
solution of this problem led to the development of parametrized Minkowski theories [6], to 
the Wigner-covariant rest-frame instant form of dynamics referred to the intrinsic inertial 
rest-frame of the N-body system [6, 7] (see Ref. [8] for Dirac's forms of the dynamics), and to 
the development of generalized radar coordinates for the synchronization of distant clocks 
in arbitrary non- inertial frames [9]. A consistent elimination of relative times (i.e. the 
choice of a convention for the synchronization of clocks, different in general from Einstein's 
convention), together with a choice of the center-of-mass time, reduce the variables of the 
N-body system to the same number as in the non-relativistic case. In special relativity 
the rest-frame instant form of dynamics, which corresponds to Einstein's synchronization 
convention due to the inertial nature of the rest-frame, implements the relativistic separation 
of the center of mass, so that only 3iV — 3 relative coordinates (or 6iV — 6 relative phase 
space variables) survive in the rest-frame. 

As it will be shown, this separation is well defined in the rest frame, which is intrinsically 
defined by the Wigner hyper-planes orthogonal to the conserved 4-momentum of the isolated 
system. It is important to stress that each such hyperplane is an instantaneous Euclidean 
and Wigner-covariant 3-space. This leads to the characterization of two distinct realizations 
of the Poincare' group, referred to an abstract observer sitting outside or inside the Wigner 
hyper-planes, namely the external and the internal realization, respectively. Also, we get the 
characterization of the relevant notions of relativistic external and internal 4- and 3- centers 
of mass, as well as of a final set of canonical relative variables with respect to the internal 
canonical 3-center of mass. In absence of interactions, such relative variables are identified 
by a canonical transformation, which, unlike the non-relativistic case, is point only in the 
momenta. However, in presence of action-at-a-distance interactions among the particles, the 
canonical transformation identifying the relative variables becomes interaction-dependent. 
In any case the use of the non-interacting internal 3-centers of mass and relative variables 
shows that the Wigner hyper-planes constitute the natural framework for the relativistic 
theory of orbits. It is also argued that a future relativistic theory of orbits will be a non- 
trivial extension of the non-relativistic theory [10], for, in the instant form of relativistic 
dynamics, the potentials appear both in the Hamiltonian and in the Lorentz boosts. 

An important point is that, as shown in Ref. [11], Jacobi normal coordinates, as well as 
notions like reduced masses and inertia tensors, do not survive in special relativity. Some 
of such notions can be recovered (and still in a non-unique way) [12] only by replacing the 
N-body system with a multi-polar expansion. A different strategy must be consequently 
devised already at the non-relativistic level [13]. 

Due to the presence of the mass-shell first-class constraints, the description of N-body 
systems in the rest-frame instant form make use of a special class of canonical transforma- 
tions, of the Shanmugadhasan type [14] and [7]. Such transformations are simultaneously 
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adapted to: i) the Dirac first-class constraints appearing in the Hamiltonian formulation 
of relativistic models (the transformations have the effect that the constraint equations are 
replaced by the vanishing of an equal number of new momenta, whose conjugate variables 
are the Abelianized gauge variables of the system); and to ii) the time-like Poincare orbits 
associated with most of their configurations. In the Darboux bases one of the final canonical 
variables is the square root of the Poincare invariant P 2 (where P M is the conserved time-like 
four- momentum of the isolated system). 

By exploiting the constructive theory of the canonical realizations of Lie groups [15, 16, 
17, 18, 19], a new family of canonical transformations was introduced in Ref.[20]. This 
family of transformations leads to the definition of the so-called canonical spin bases, in 
which also the Pauli-Lubanski Poincare invariant W 2 = —P 2 S T for the time-like Poincare 
orbits becomes one of the final canonical variables (provided the rest-frame Thomas spin 
St = Yl a =i S a is different from zero). The essential point is that the construction of the 
spin bases exploits the clustering of spins instead of the Jacobi clustering of centers of mass, 
which is an ill defined notion at the relativistic level. 

Note that, in spite of its genesis in a relativistic context, the technique used in the 
determination of the spin bases, related to a typical form [15] of the canonical realizations 
of the E(3) group, can be easily adapted to the non-relativistic case, where W 2 is replaced 
by the invariant S 2 of the extended Galilei group. The fact that the traditional Jacobi 
clustering of the centers of mass of the sub-clusters is replaced by the clustering of the spins 
S a , (a — 1, .., N— 1) of the sub-clusters, as in the composition of quantum mechanical angular 
momenta, is the basic trick that makes the treatment of the non-relativistic N-body problem 
directly extendible to the relativistic case. The clustering can be achieved by means of suitable 
canonical transformations, which in general are non-point both in the coordinates and the 
momenta. This entails that the quantum implementations of such canonical transformations 
as unitary transformations be non-trivial. The extension of our formalism to quantum 
mechanics remains an open problem. 

Our aim at this point the construction of a canonical Darboux basis adapted to the non- 
Abelian SO (3) symmetry, both at the non-relativistic and the relativistic level. The three 
non-Abelian Noether constants of motion S = S r f r (f r are the axes of the inertial laboratory 
or space frame) are arranged in these canonical Darboux bases as an array containing the 
canonical pair S 3 , (3 = tg' 1 ^ and the unpaired variable S — \S\ (scheme A of the canonical 

realization of SO(3) [16]; the configurations with S = are singular and have to be treated 
separately). The angle canonically conjugated to S, say a, is an orientational variable, 
which, being coupled to the internal shape degrees of freedom, cannot be a constant of 
motion. However, being conjugated to a constant of the motion, it is an ignorable variable 
in the Hamiltonian formalism, so that its equation of motion can be solved by quadratures 
after the solution of the other equations. In conclusion, in this non-Abelian case one has 
only two (instead of three as in the Abelian case) commuting constants of motion, namely 
S and S* 3 (like in quantum mechanics). This is also the outcome of the momentum map 
canonical reduction [21, 22] by means of adapted coordinates. Let us stress that a, S 3 , (3 
are a local coordinatization of any co- adjoint orbit of SO (3) contained in the N-body phase 
space. Each co-adjoint orbit is a 3-dimensional embedded sub- manifold and is endowed with 
a Poisson structure whose neutral element is a. By fixing non-zero values of the variables S 3 , 
(3 = tg' 1 ^ through second-class constraints, one can define a (6N-8)-dimensional reduced 
phase space. The canonical reduction cannot be furthered by eliminating S, just because a is 
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not a constant of motion. Yet, the angle a allows us to construct a unit vector R, orthogonal 
to S, such that S, R, S x R (the notation "means unit vector) is an orthonormal frame that 
we call spin frame. 

The final lacking ingredient for our construction of body frames comes from the group- 
theoretical treatment of rigid bodies [22] (Chapter IV, Section 10). Such treatment is based 
on the existence of realizations of the (free and transitive) left and right Hamiltonian actions 
of the SO (3) rotation group on either the tangent or cotangent bundle over their configura- 
tion space. The generators of the left Hamiltonian action [23], which is a symmetry action, 
are the above non-Abelian constants of motion S 1 , S 2 , S 3 , [{S r , S s } = e rsu S u ]. 

At this point let us stress that, in the approach of Ref. [1] the SO(3) principal bundle is 
built starting from the relative configuration space and, upon the choice of a body-frame 
convention, a gauge- dependent SO (3) right action is introduced. The corresponding task 
in our case is the following: taking into account the relative phase space of any isolated 
system, we have to find out whether one or more SO (3) right Hamiltonian actions could be 
implemented besides the global SO (3) left Hamiltonian action. In other words, we have to 
look for solutions S r , r=l,2,3, [with Er(^ r ) 2 = Er(^ r ) 2 = S 2 ], o f the partial differential 
equations {S r , S s } = 0, {S r , S s } = — e rsu S u and then build corresponding left invariant 
Hamiltonian vector fields. Alternatively, one may search for the existence of a pair S 3 , 
7 = of canonical variables satisfying {7, S* 3 } = —1, {7, S r } = {S 3 ,S r } = and 

also {7, a} = {S 3 ,a} = 0. Local theorems given in Refs.[15, 16] guarantee that this is 
always possible provided A > 3. Clearly, the functions S r , which are not constants of the 
motion, do not generate symmetry actions. What matters here is that each explicitly given 
right action leads to the characterization of the following two structures: i) a dynamical 
reference frame (say N, x, N x x), that we call dynamical body frame; ii) a canonical spin 
basis including both the orientational and the shape variables. 

In conclusion, we show that, after the center-of-mass separation, by exploiting the new 
notions of dynamical body frames and canonical spin bases, it is possible to build a geo- 
metrical and group-theoretical procedure for the common characterization of the rotational 
kinematics of non-relativistic and relativistic N-body systems. The two cases are treated in 
Sections II and III, respectively. 

In the last Subsection of Section III we show that the relativistic separation of the center 
of mass, realized by the rest-frame instant form, can be extended to continuous deformable 
relativistic isolated systems, namely relativistic field configurations, strings and fluids. The 
action principle of such systems can be transformed into a parametrized Minkowski theory on 
space-like hyper-surfaces whose embeddings in Minkowski space-time are the gauge variables 
connected with the arbitrariness in the choice of the convention for clock synchronization 
(namely the choice of the equal-time Cauchy surfaces for the field equations). Then the rest- 
frame instant form on Wigner hyper-planes emerges in a natural way also for fields (ADM 
canonical metric [24] and tetrad gravity [7] naturally deparametrize to it). All the notions 
of external and internal 4- and 3-centers of mass can be extended to field configurations 
under the condition that a collective 4- vector [25], canonically conjugate to the configuration 
conserved 4-momentum, be definable. 

The construction of such a collective variable, with respect to which the energy- 
momentum distribution of the configuration itself is peaked, is exemplified for a classical real 
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Klein-Gordon field [26]. New features, like an internal time variable, absent in the particle 
case, emerge for fields and entail that each constant energy surface of the configuration be a 
disjoint union of symplectic manifolds. These results can be extended to relativistic perfect 
fluids [27]. 

The next issue has to do with the simulation of an extended system by means of as 
few as possible global parameters, sufficient to maintain an acceptable phenomenological 
description of the system. This can be achieved by replacing the extended system with a 
multipolar expansion around a world-line describing its mean motion. After many attempts, 
a general approach to this problem has been given by Dixon in Ref. [28] for special relativity 
and in Ref. [29] (see also Refs.[30, 31]) for general relativity, after a treatment of the non- 
relativistic case [29]. 

With this in view, we show in Section IV that the rest-frame instant-form of the dynamics 
is the natural framework (extendible to general relativity [7]) for the formulation of rela- 
tivistic multipolar expansions [12] exhibiting a clear identification of the internal canonical 
3-center of mass (with associated spin dipole and higher multi-poles) as a preferred center of 
motion. Besides a system of N free relativistic particles, we also discuss an open system de- 
fined by cluster of n < N charged particles inside an isolated system of N charged particles, 
plus the electro-magnetic field in the radiation gauge [6]. Finally, as a prerequisite to the 
treatment of relativistic perfect fluids [27], we give some results concerning a configuration of 
a classical Klein-Gordon field [26], where the collective variable allows identifying a natural 
center of motion together with the associated Dixon multi-poles. 

II. THE NON-RELATIVISTIC CANONICAL SPIN BASES AND DYNAMICAL 
BODY FRAMES. 

Let us consider N free non- relativistic particles of masses m^, % = 1,..,N, described by 
the configuration variables ffi and by the momenta «j. The Hamiltonian H = ^2f =1 is 

restricted by the three first-class constraints k + = J^li ^ ~ defining the center-of-mass 
or rest frame. 

Let us introduce the following family of point canonical transformations (m = YliLi m «) 
realizing the separation of the center of mass from arbitrary relative variables (for instance 
some set of Jacobi normal coordinates) 

—* 

defined by: 

1 Ar_1 m N ~ l 

* / AT m 





Pa 


«+ 





(2.1) 
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N N 

ST^ THi _ _ _ \ ^ _ 

m 

t=l i=l 



N I N iV_1 

p a = V^^'JaiVi, K a = —= \2 T ai Ki, S = Y] p a X 7T a 

1=1 v i=l a=l 



AT „ N 



r - -r -lw 

J- ai Tai / . Tafc) Taj J- ai at / , afc> 

k=l k=l 
N N 

Elai = 0, V — T ai = 0, 



. m 

i=i i=i 

TV AT 

i=l i=l 

JV-1 j AT-1 

7ai7aj = Sij -—, Tai^ = 5ij -. (2.2) 

iV 771 
a=l a=l 

Here, the 7 ai 's and the r ai 's are numerical parameters depending on |(iV — l)(iV — 2) free 
parameters [6]. 

It can be shown [13] that the relative motion is described by the following Lagrangian 
and Hamiltonian 



j 1..JV-1 

A-ej(f) = - Y k ab[mi,T ai \p a (t) ■ p b (t) , 



2 

a,b 



k a b[ m i,T C i] — ki, a [mi,T C i\ — — ^ ^ mJaJfo, 



i=l 

AT 



,=1 m * 



4 

AT-1 

^<»(*) = ^2 k ab[mi,r c i\Pb(t), 



b=l 

1..N-1 



H rel = \ ^2 k ab[ m i^ r ai\ K a (t) •??&(*)■ (2-3) 



2 

ab 



If we add the gauge fixings ~ and we go to Dirac brackets with respect to the second- 
class constraints k + ~ 0, ~ 0, we get a (6N-6)-dimensional reduced phase space spanned 
by p a , and with S = J2 a =i Pa x n a = J2 a =i S a - At the non-relativistic level [1] the next 
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problem of the standard approach for each N is the diagonalization of the matrix k a b[mi, r o j]. 
The off-diagonal terms of the matrix k a b[mi,T a i\ are called mass polarization terms, while 
its eigenvalues are the reduced masses (see for instance Ref. [32]). In this way the Jacobi 
normal coordinates p a = s a , with conjugate momenta n a = n sa , are introduced. 

In the rest frame the 11 generators of the extended Galilei group (the total mass m is a 
central charge) [17] are 

N N ^ 9 N 

K7 



1=1 1=1 1=1 



m 

N N-l N-l 



J — ^2 rfi x ^i- S S a = ^2 p a x n a , 

i=l a=l a=l 

N 

K = -^2 m iVi + K+t w -mq nr , (2.4) 



i=l 



and the gauge fixings q nr ~ are equivalent to K fa 0. 

Following the strategy delineated in the Introduction, we search for a geometrical and 
group-theoretical characterization of further canonical transformations leading to a privi- 
leged class of canonical Darboux bases for the N-body. Specifically, they must be adapted 
to the SO(3) subgroup [15, 16] of the extended Galilei group and contain one of its invariants, 
namely the modulus of the spin: 

1) Every such basis must be a scheme B (i.e. a canonical completion of scheme A, 
according to the language of Ref. [15, 16, 17, 20]) for the canonical realization of the rotation 
group SO(3), viz. it must contain its invariant S and the canonical pair S 3 , (3 = tg -1 ^. 
This entails that, except for a, all the remaining variables in the canonical basis be SO(3) 
scalars. 

2) The existence of the angle a satisfying {a, S} = 1 and {a, S 3 } = {a, f3} = leads 

^ — * 

to the geometrical identification of a unit vector R orthogonal to S and, therefore, of an 
orthonormal frame, the spin frame S, R, S x R. 

3) From the equations R 2 = 1 and {S ■ R,R 1 } = it follows the symplectic result 
{R\ R?} = 0. As a byproduct, we get a canonical realization of a Euclidean group E(3) with 
generators S, R [{R\ R j } = 0, {R\ S j } = e ijk R k ] and fixed values of its invariants R 2 = 1, 
S ■ R = (non-irreducible type 3 realization according to Ref. [20]). 

4) In order to implement a SO (3) Hamiltonian right action in analogy with the rigid 
body theory [22], we must construct an orthonormal triad or body frame N, x, iV x \- The 
decomposition 

S = S 1 X + S 2 N x X + S 3 N = f S r e r , (2.5) 

identifies the SO(3) scalar generators S r of the right action provided they satisfy {S r , S s } = 
—e rsu S u . This latter condition together with the obvious requirement that N, x, N x x be 
SO(3) vectors [{iV r ,5 s } = e rsu N u , {x r ,S s } = e rsu X u , {N x x , S s } = e rsu N x £] entails 
the equations 
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{A> N s } = {A> x s } = {X r , X s } = 0. (2.6) 

Each solution of these equations identifies a couple of canonical realizations of the E(3) 
group (non-irreducible, type 2): one with generators S, N and non- fixed invariants S 3 = S-N 
and \N\; another with generators S, x an d non-fixed invariants S 1 = S ■ x an d \x\- Such 
realizations contain the relevant information for constructing the new canonical pair S 3 , 
7 = tg~ 1 ^: of SO(3) scalars. Since {a, S 3 } = {0,7} = must hold, it follows [20] that 
the vector N necessarily belongs to the S-R plane. The three canonical pairs S, a, S 3 , (3, 
S 3 , 7 will describe the orientational variables of our Darboux basis, while \N\ and \x\ will 
belong to the shape variables. For each independent right action [i.e. for each solution N, 
X of Eqs.(2.6)], we can identify a canonical spin basis containing the above 6 orientational 
variables and 6iV — 12 canonical shape variables. Alternatively, an anholonomic basis can be 
constructed by replacing the above six variables by S r (or S r ) and three uniquely determined 
Euler angles a, (3, 7 (see Ref. [13]). Let us stress that the non-conservation of S r entails that 
the dynamical body frame evolves in a way dictated by the equations of motion, just as it 
happens in the rigid body case. 

We can conclude that the N-body problem has hidden structures leading to the charac- 
terization of special dynamical body frames which, being independent of gauge conditions, 
are endowed with a direct physical meaning. 

A. The 2-body system. 

For N = 2, a single E(3) group can be defined: it allows the construction of an orthonor- 
mal spin frame S, R, R x S in terms of the measurable relative coordinates and momenta 

.— — -»2 

of the particles. The relative variables are p = p , 7f and the Hamiltonian is H re \ = |-, 

where p = ™^™^ 2 is the reduced mass. The spin is S = p x n. Let us define the following 
decomposition 

— # 

p = pR, P = \[f, R = -=P, R 2 = l, 

P 

S - S 

7T = TlR R X S = Tip f) X S, 

P ' P 

— * 

n = jf.R = jf.p, £ = 5-^ = 0. (2.7) 

As shown in Ref. [20], it is instrumental considering the following non-point canonical 

— # 

transformation adapted to the SO (3) group, valid when S ^ 



a = tg-^(p-n-^n 3 ). (2.8) 

In the language of Ref. [20], the two pairs of canonical variables a, S, (3, S 3 form the 
irreducible kernel of the scheme A of a (non-irreducible, type 3, ) canonical realization of 



P 



a (3 


P 


S S 3 


TT 
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the group E(3), generated by S, R, with fixed values of the invariants R 2 = 1, R ■ S = 0, 
just as the variables S 3 , (3 and S form the scheme A of the SO(3) group with invariant S. 
Geometrically, we have: i) the angle a is the angle between the plane determined by S and 
/ 3 and the plane determined by S and R; ii) the angle (3 is the angle between the plane S 
- f 3 and the plane / 3 - f ± . Moreover, S 1 = ^/(S) 2 - (S 3 ) 2 cos (3, S 2 = ^/(S) 2 - {S 3 ) 2 sin(3, 
R 1 = p 1 = sinf3sina — ^-cos (3cos a, R 2 = p 2 = —cosfisina — ^-sin (3cos a, R 3 = p 3 = 

In this degenerate case (N=2), the dynamical shape variables p, tc coincide with the static 
ones and describe the vibration of the dipole. The rest-frame Hamiltonian for the relative 

motion becomes (/ is the barycentric tensor of inertia of the dipole) H re i = \ I~ 1 S 2 + , 
/ = pp 2 , while the body frame angular velocity is u - (Ul "' - - 



dS I ■ 



B. The 3-body system. 

—*—*—» — » 
For N — 3, where we have S = Si + S 2 , a pair of E(3) groups emerge, associated with Si 

— * A, 

and respectively. We have now two unit vectors R a and two E(3) realizations generated 
by S a , R a respectively and fixed invariants R 2 — 1, S a • R a — (non-irreducible, type 2, 
see Ref. [20]). We shall assume S ^ 0, because the exceptional SO(3) orbit S = has to be 
studied separately by adding S ~ as a first-class constraint. 

For each value of a = 1, 2, we consider the non-point canonical transformation (2.8) 



Pa 



da Pa 


Pa 


S a S 3 


Via 



Pa — PaRa, Pa — Ra — — — Pa-, — 1, 

Pa 

s ~ 

7T a = TX a R a + — S a X ,R a , 7T a = 7T a • i? a . (2.9) 

Pa 

In this case, besides the orthonormal spin frame, an orthonormal dynamical body frame 
N, x, N x x, i.e. a SO (3) Hamiltonian ng/zt action, can be defined. The simplest choice, 
within the existing arbitrariness [34], for the orthonormal vectors iV and x functions only of 
the relative coordinates is 



N = \ { Ri + R 2 ) = \ { pi + p 2) , N=^. |rf| = -' 1 + *-* 



x = \(Ri-R2) = \(pi-P2), x=||, 1x1 = / i l ' p2 = ^Ji-N 2 , 
i , - . l 



Nxx = -2P1XP2, |iVxx| = |iV||x| = 2Vl-(Pi-P2) 2 , N-x = 0. (2.10) 
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As said above, (Eq.(2.5), this choice is equivalent to the determination of the non- 
conserved generators S r of a Hamiltonian right action of SO (3). 

The realization of the E(3) group with generators S, A and non-fixed invariants A 2 , S-N 
leads to the final canonical transformation introduced in Ref. [20] 



Pa 



Oil Pi 


Ol2 132 


Pa 


Si Sf 


S2 S 3 





a (3 7 


\N\ 


Pat 


S = S S 3 S 3 = S- 


A 







[1 - Pi ■ P2W 1 + Pi' P2 



(2.11) 



For N=3 the dynamical shape variables, functions of the relative coordinates p a only, are 
\N\ and p a , while the conjugate shape momenta are £, n a . 

— * /\ — * 

We can now reconstruct S and define a new unit vector R orthogonal to S by adopting the 
prescription given after Eq.(2.8). The vectors S , R, S x R build up the spin frame for N=3. 



The angle a conjugate to S is explicitly given by a = —tg~ 



(SxN) s 



-tg- 



(SxRf 



[Sx(SxN)] 3 " a [Sx(SxR)] 3 ' 

The two expressions of a given here are consistent with the fact that S, R and A are 
coplanar, so that R and N differ only by a term in S. 

As a consequence of this definition ofR, we get the following expressions for the dynamical 
body frame N, x, N x x m terms of the final canonical variables 



A 



S 3 



1 



X 



cosi)S + sini>R= — S + -J(S) 2 - (S 3 ) 2 R = 
A[S,a;S 3 ,/3;S 3 , 7 ], 

sin ipcos '-/S — cos ipcos ^R + sin x R — 

o ~ S 3 S^ R + S 2 S x R o -o , 

~s s - -j / =- = xfr 5 > # 5 > 7], 



(S) 2 - (S 3 ) 2 



4 



5 



R 



Rx S 



= sin ipcos 7X + sin ipsin 7 A x x + cos 

^ 1 [5^ + 5 2 A x x + S 3 A , 

= — cos ipcos 'yx ~ cos ipsin x x + sinipN, 

= — sin 7X + cos 7 A x x- 



(2.12) 



While ip is the angle between S and A, 7 is the angle between the plane A — x and the 
plane S — A. As in the case N=2, a is the angle between the plane S — K and the plane 
S — R, while (3 is the angle between the plane S — f 3 and the plane f 3 — fi- 
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Finally, as shown in Ref. [13], we can perform a sequence of a canonical transformation 
to Euler angles a, (3, 7 (explicitly calculable in terms of p a and 7T a ) with their conjugate 
momenta, followed by a transition to the anholonomic basis used in the orient at ion- shape 
bundle approach [1] [cf = (pi,p2, |^|), Pn = (^1,^2, are the dynamical shape variables) 



a (3 7 


\N\ 


Pa 


S = S S 3 S 3 


t 


Via 



a (3 7 


|iV| 


Pa 


S 1 S 2 S 3 







a /3 7 


q^(Pa) 


s 1 s 2 s 3 


Pli,(Pa, TTa) 



(2.13) 



The Euler angles a, (3, 7 are determined as the unique set of dynamical orientation 
variables. 

For JV = 3 we get a result similar to the orient at ion- shape bundle approach, namely 
the following relation between the space and body components of the relative coordinates: 

p' a = KU&,M)f? a (q), with pl(q) = (-) a+1 p a Vl - N\ p 2 a (q) = 0, p 3 a (q) = p a \N\, S r = 
lZ r s (a, (3, 7)S S . It can also be shown that for N=3 our definition of dynamical body frame 
can be reinterpreted as a special global cross section (xxzz gauge, where x stays for x an d 
z for N; an outcome that is independent of the particular choice made for TV and x) °f the 
trivial SO (3) principal bundle of Ref. [1] , namely a privileged choice of body frame. While the 
above dynamical body frame can be identified with the global cross section corresponding to 
the xxzz gauge, all other global cross sections cannot be interpreted as dynamical body frames 
(or dynamical right actions), because the SO(3) principal bundle of Ref. [1] is built starting 
from the relative configuration space and, therefore, it is a static, velocity-independent, 
construction. As a matter of fact, after the choice of the shape configuration variables g M 
and of a space frame in which the relative variables have components p r a , the approach 
of Ref. [1] begins with the definition of the body-frame components p r a {q^) of the relative 
coordinates, in the form p r a = R rs (8 a )p s a (q' J '), and then extends it in a velocity-independent 

def 

way to the relative velocities p r a = R rs (6 a )v s a . Here R is a rotation matrix, 6 a are arbitrary 
gauge orientational parameters and p r a {q >J ') is assumed to depend on the shape variables only. 
In our construction we get instead p r a = R rs (a, /3, 7)p a ((? M ) in the xxzz gauge, so that in 
the present case (N=3) all dynamical variables of our construction coincide with the static 
variables in the xxzz gauge. On the other hand, in the relative phase space, the construction 
of the evolving dynamical body frame is based on non-point canonical transformation. 



C. The N-body system. 

For N=4, it holds S = Si + S 2 + S3. Since we have three relative variables pi, p 2 , f>3 
and related momenta 7?i, 712, vr3, it is possible to construct three sets of spin frames and 
dynamical body frames corresponding to the hierarchy of clusterings ((ab)c) [i.e. ((12)3), 
((23)1), ((31)2); the subscripts a, b, c denote any permutation of 1, 2, 3] of the relative spins 
S a . The associated three canonical Darboux bases share the three variables S 3 , (3, S (viz. 
S), while both the remaining three orientational variables and the shape variables depend 
on the spin clustering. This entails the existence of three different SO (3) right actions 
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with non-conserved canonical generators S? A y A=l,2,3. Consistently, one can define three 

anholonomic bases <S(a), 13(A), 7(A), Sfy and associated shape variables q*? A y P(A)n, P — 1, •-, 6, 
connected by canonical transformations leaving S r fixed. The relative variables are therefore 
naturally split in three different ways into 6 dynamical rotational variables and 12 generalized 
dynamical shape variables. Consequently, we get three possible definitions of dynamical 
vibrations. 

By using the explicit construction given in Appendix C of Ref. [13], we define the following 
sequence of canonical transformations (we assume S ^ 0; Sa ^ 0, A = a, b, c) corresponding 
to the spin clustering pattern abc t— > (ab)c t— > ((ab)c) [build first the spin cluster (ab), then 
the spin cluster ((ab)c)]: 



Pa Pb Pc 
VT a ^b 









Pa Pb Pc 




s b s 3 


s c s 3 c 


K b TT C 



(ab)c 



tt(ab) P(ab) 7(ab) 




=^ 

N( ab )| Pa Pb Pc 




^(ab) ^fafe) ^fab) = ^at) ' N {«*>) 




£(ab) K a lT b TT C 


Ot((ab)c) P({ab)c) 1((ab)c) 


-i =? 

A^((a6) c ) 7(o6) 


=? 

\N( ab ) p a p b Pc 


s = s s 3 s 3 = s ■ % 6)C ) 


£((ab)c) S'(ab) • A r (afe) £(ab) ^6 7T C 



non can. 



a (3 7 


— * — * 

A"(( ab ) c ) 7( a6 ) ^(ab) 


Pa Pb Pc 


S l S 2 S 3 


£((ab)c) fya6) = <S(a&) • N( ab ) £( o6 ) 7T a 7T 6 7T C 



(2.14) 

The first non-point canonical transformation is based on the existence of the three unit 
vectors Ra, A = a,b, c, and of three E(3) realizations with generators Sa, Ra and fixed 
values (R 2 A — 1, Sa • Ra — 0) of the invariants. 

In the next canonical transformation the spins of the relative particles a and b are coupled 
to form the spin cluster (ab), leaving the relative particle c as a spectator. We use Eq.(2.10) 
to define N {ab) = \(R a + Rb), X(ab) = \(R a ~ Rb), S {ab) = S a + S b , W {ab) = S a - S b . We get 
N {ab) ■ X(ab) = 0, {N{ ab) , N° ab) } = {N r {ab) , X s {ab) } = {x\ ab) , X s {ab) } = and a new E(3) realization 

generated by S( ab ) and N^, with non-fixed invariants \N( ab )\, S( ab ) ■ N( ab ) = fi( a b)- It can be 
shown that it holds 



Pa = Pa 



Pb = Pb 



\N {a b)\N iab) + 



N 2 



(ab)^(ab) 



\N(ab)\N(ab) 



N 2 



(aft) X (aft) 



Pc — PcRc- 



(2.15) 



It is then possible to define a^ ab ) and P( ab ) and a unit vector R^ ab ) with 5( ob ) ■ R( a b) = 0, 
{R[ a b), R S (ab)} = 0- This unit vector identifies the spin cluster (ab) in the same way as the 
unit vectors Ra = Pa identify the relative particles A. 
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The next step is the coupling of the spin cluster (ab) with unit vector R( a b) [described by 
the canonical variables a^ ab ), S( ab ), (3( ab ) Sfab)] ^° ^ ne re ^ a ^ ve particle c with unit vector R c 
and described by a c , S c , f3 c , S^: this characterizes the spin cluster ((ab)c). 

Again, Eq.(2.10) allows to define N {{ab)c) = \(R(ab) + Rc), X((ab)c) = \(R{ab) ~ Rc), 
S = S(( ab)c) = S(ab) + S c , W^ ab)c) = S( ab ) - S c . Since we have N i{ab)c) ■ X{(ab)c) = and 

{ N ((ab)c)i N ((ab)c)} = { N ((ab)c)iX((ab)c)} = {X((a6)c) ' X(( ab ) c ) } = due to {R r {ab) , R*} = 0, 

a new E(3) realization generated by S and N^ ab ) c ) with non-fixed invariants |iV(( a j,) c )|, 
S ■ N((ab)c) = "S* 3 emerges. Then, we can define a^ ab ) c ) and /3(( &) c ) an d identify a final 
unit vector ^ ((ab)c) with S ■ R {{ab)c) = and {R r ({ab)c) , R s {{ab)c) } = 0. 

In conclusion, when S ^ 0, we find both a spin frame S, R(( a b)c), R{(ab)c) x S and a 
dynamical body frame X{(ab)c), N ({ab)c) X X{(ab)c), N^ ab)c) , like in the 3-body case. There is an 
important difference, however: the orthonormal vectors N(( a b)c) and X({ab)c) depend on the 
momenta of the relative particles a and b through R( a b), so that our results do not share any 
relation with the N=4 non-trivial SO (3) principal bundle of the orientation-shape bundle 
approach. 

The final 6 dynamical shape variables are g M = {\N^ ab ^\,^f^, \N^\, p a , p b , p c }. While 
the last four depend only on the original relative coordinates pa, A — a,b, c, the first two 
depend also on the original momenta tta'- therefore they are generalized shape variables. In 
Ref.[13] it is shown that, instead of p r a = R rs (a, f3, l)p s a (q^), it holds 

p r A = TZ rs (a, 0, 7) p s A (q", p„ S r ), A = a, b, c. (2.16) 

Clearly, this result stands completely outside the orientation-shape bundle approach. As 
a consequence the above anholonomic bases and the associated evolving dynamical body 
frames, however, have no relations with the N=4 static non-trivial SO (3) principal bundle 
of Ref.[l], which admits only local cross sections. Each set of 12 generalized dynamical 
canonical shape variables is obviously defined modulo canonical transformations so that it 
should even be possible to find local canonical bases corresponding to the local cross sections 
of the N=4 static non-trivial SO(3) principal bundle of Ref.[l]. 

Finally, it can be shown that, starting from the Hamiltonian H re i^ ab ^ expressed in the fi- 
nal variables, we can define a rotational Hamiltonian H^) , and a vibrational Hamiltonian 

' rel({ab)c) 

H rei{{ab)c) (vanishing of the physical dynamical angular velocity ^'{ {ab)c) = 0), but i? re /((o6) c ) 
fails to be the sum of these two Hamiltonians showing once again the non-separability of 
rotations and vibrations. Let us stress that in the rotational Hamiltonian we find an inertia- 
like tensor depending only on the dynamical shape variables. A similar result, however, does 
not hold for the spin-angular velocity relation. 

The price to be paid for the existence of 3 global dynamical body frames for N=4 is a 
more complicated form of the Hamiltonian kinetic energy. On the other hand, dynamical 
vibrations and dynamical angular velocity are measurable quantities in each dynamical body 
frame. 

Our results can be extended to arbitrary N, with S = J2^=i ^a- There are as many 
independent ways (say K) of spin clustering patterns as in quantum mechanics. For instance 
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for N=5, K — 15 : 12 spin clusterings correspond to the pattern (((ab)c)d) and 3 to the 
pattern ((ab)(cd)) [a,b,c,d = 1, ..,4]. For N=6, K = 105: 60 spin clusterings correspond to 
the pattern ((((ab)c)d)e), 15 to the pattern (((ab)(cd)e) and 30 to the pattern (((ab)c)(de)) 
[a,b,c,d,e = 1,..,5]. Each spin clustering is associated to: a) a related spin frame; b) a 

related dynamical body frame; c) a related Darboux spin canonical basis with orientational 

s 2 

variables S 3 , [3, S, ot(A), S?a)i 7(A) — ^g~ x jr^-, A — 1, ..,K (their anholonomic counterparts 

are (5(a), P(a), 7(A), S(a.) w ith uniquely determined orientation angles) and shape variables 
5(A)' £V(A)> P — 1, ■■■> 3iV — 6. Furthermore, for N > 4 we find the following relation between 
spin and angular velocity: S r = 1 rs {q[ A) ) u[ A) + /^(^(A)^ 

Therefore, for iV > 4 our sequence of canonical and non-canonical transformations leads 
to the following result, to be compared with Eq.(2.13) of the 3-body case 



Pa 



tta 



a (3 7 


g M (PA,VTA) 


s 1 s 2 s 3 


P^Pa, v?a) 



(2.17) 



Therefore, for N > 4 and with S ^ 0, Sa ^ 0, A = a, b, c, namely when the standard 
(3N-3)-dimensional orientation-shape bundle is not trivial, the original (6N-6)-dimensional 
relative phase space admits as many dynamical body frames as spin canonical bases, which 
are globally defined (apart isolated coordinate singularities) for the non-singular N-body con- 
figurations with S 7^ (and with non-zero spin for each spin sub-cluster). These dynamical 
body frames do not correspond to local cross sections of the static non-trivial orientation- 
shape SO (3) principal bundle and the spin canonical bases do not coincide with the canonical 
bases associated to the static theory. 

— * 

The 5 = 0, C-horizontal, cross section of the static SO (3) principal bundle corresponds 
to N-body configurations that cannot be included in the previous Hamiltonian construction 
based on the canonical realizations of SO (3): these configurations (which include the singular 
ones) have to be analyzed independently since they are related to the exceptional orbit of 
SO (3), whose little group is the whole group. 

While physical observables must be obviously independent of the gauge- dependent static 
body frames, they do depend on the dynamical body frame, whose axes are operationally 
defined in terms of the relative coordinates and momenta of the particles. In particular, a 
dynamical definition of vibration, which replaces the S = 0, C-horizontal, cross section of the 
static approach [1], is based on the requirement that the components of the angular velocity 
vanish. Actually, the angular velocities with respect to the dynamical body frames become 
now measurable quantities, in agreement with the phenomenology of extended deformable 
bodies (see, e.g., the treatment of spinning stars in astrophysics). 

In this connection, let us recall that the main efforts done in developing canonical trans- 
formations for the N-body problem have been done in the context of celestial mechanics. As 
a consequence, these transformations are necessarily adapted to the Newtonian gravitational 
potential. As well known, the Kepler problem and the harmonic potential are the only inter- 
actions admitting extra dynamical symmetries besides the rotational one. The particularity 
of the final canonical transformations which have been worked out in that mechanical con- 
text is that the Hamiltonian figures as one of the final momenta. This is in particular true 
for the N=2 body problem (Delaunay variables; our f3 is the longitude of the ascending node 
Q [10]) and for the N=3 body problem (Jacobi's method of the elimination of the nodes). 
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On the other hand, it is well-known that, for generic non-integrable interactions, putting 
the Hamiltonian in the canonical bases is quite useless, since it does not bring to isolating 
integrals of the motion having the capability of reducing the dimensionality of phase space. 
Here we stress again that our construction rests only on the left and right actions of the 
SO(3) group and is therefore completely independent of form of the interactions. 

The expression of the relative Hamiltonian H re i for N = 3 in these variables is a compli- 
cated function, given in Ref. [13], which fails to be the sum of a rotational plus a vibrational 
part. It is an open problem for N=3 and iV > 4 to check whether suitable choices of the 
numerical constants 7 ai , more general body frames obtained by exploiting the freedom of 
making arbitrary configuration-dependent rotations [34] and/or suitable canonical transfor- 
mations of the shape variables among themselves can simplify the free Hamiltonian and/or 
some type of interaction. 

In conclusion, for each N, we get a finite number of physically well-defined separations 
between rotational and vibrational degrees of freedom. The unique body frame of rigid 
bodies is replaced here by a discrete number of evolving dynamical body frames and of spin 
canonical bases. Both of them are grounded on patterns of spin couplings which are the 
direct analogues of the coupling of quantum-mechanical angular momenta. 

III. THE RELATIVISTIC CENTER-OF-MASS PROBLEM, RELATIVE VARI- 
ABLES AND THE ROTATIONAL KINEMATICS IN SPECIAL RELATIVITY. 

Let us see what happens when we replace Galilean space plus time with Minkowski 
space-time, already for the simple model-system of N free scalar positive-energy particles. 

First of all we have to describe a relativistic scalar particle. Among the various possi- 
bilities (see Refs.[6, 35] for a review of the various options) we will choose the manifestly 
Lorentz covariant approach based on Dirac's first-class constraints 

p^-em^O. (3.1) 

The associated Lagrangian description is based on the 4- vector positions xf(r) and the 

action S = f dr(^ — e £V rai^J exj{r) j , where r is a Lorentz scalar mathematical time, i.e. 

an afline parameter for the particle time-like world-lines. Then, Lorentz covariance implies 
singular Lagrangians and the associated Dirac's theory of constraints for the Hamiltonian 
description. The individual time variables x°(t) are the gauge variables associated to the 
mass-shell constraints, which have the two topologically disjoint solutions p° ~ ±^/mf + . 
As discussed in Ref. [36] and [6] this implies that: 

i) a combination of the time variables can be identified with the clock of one arbitrary 
observer labeling the evolution of the isolated system; 

ii) the N — 1 relative times are related to observer's freedom of looking at the N particles 
either at the same time or with any prescribed relative delay, or, in other words, to the 
convention for synchronization of distant clocks (definition of equal-time Cauchy surfaces on 
which particle's clocks are synchronized) used by the observer to characterize the temporal 
evolution of the particles [9]. 

Introducing interactions in this picture without destroying the first-class nature of the 
constraints is a well-known difficult problem, reviewed in Refs.[6, 35], where also the models 
with second-class constraints are considered and compared. 



18 



A. Parametrized Minkowski theories. 

If the particle is charged and interacts with a dynamical electromagnetic field, a problem 
of covariance appears. The standard description is based on the action 

S = -em J dr^ex 2 (T)-e J dr J d 4 z5 4 {z - x^x^A^z) - ^ J d 4 zF^(z)F^(z). 

(3.2) 

By evaluating the canonical momenta of the isolated system, charged particle plus electro- 
magnetic field, we find two primary constraints: 

X(r) = {p-eA(x(T)) S j 2 -em 2 « 0, n°(z°, z) « 0. (3.3) 

It is immediately seen that, since there is no concept of equal time, it is impossible to 
evaluate the Poisson bracket of these constraints. Also, due to the same reason, the Dirac 
Hamiltonian, which would be H D = H c + A(t)x(t) + / d 3 z\°(z°, z)n (z , z) with H c the 
canonical Hamiltonian and with A(r), X°(z°, z) Dirac's multipliers, does not make sense. 
This problem is present even at the level of the Euler-Lagrange equations, specifically in 
the formulation of a Cauchy problem for a system of coupled equations some of which are 
ordinary differential equations in the afline parameter r along the particle world-line, while 
the others are partial differential equations depending on Minkowski coordinates z^. Since 
the problem is due the lack of a covariant concept of equal time between field and particle 
variables, a new formulation is needed. 

In Ref.[6], after a discussion of the many-time formalism, a solution of the problem was 
found within a context suited to incorporate the gravitational field. The starting point is an 
arbitrary 3+1 splitting of Minkowski space-time with space-like hyper-surfaces (see Ref.[9] 
for more details on the admissible splittings). After choosing the world-line of an arbitrary 
observer, a centroid x%(t), as origin, a set of generalized radar 4-coordinates [9], adapted to 
the splitting, is provided by the r afline parameter of the centroid world-line together with 

a system of curvilinear 3-coordinates a = {a r } vanishing on the world-line. In this way 
we get an arbitrary extended non-inertial frame centered on the (in general accelerated) 
observer described by the centroid. The coordinates (r, a) are generalized radar coordinates 
depending upon the choice of the centroid and the splitting. The space-like hyper-surfaces 
are described by their embedding z^(r,a) in Minkowski space-time. The metric induced 
by the change of coordinates x^ i— > a A = (r, a) is <7ab(t, <?) = 8a z^(r, a) r/^ ds z u (r,a) 
(Oaz^It, a) = d z^(r,a) / do~ A ). This is essentially Dirac's reformulation [37] of classical field 
theory (suitably extended to particles) on arbitrary space-like hyper- surfaces {equal time or 
simultaneity Cauchy surfaces). Note, incidentally, that it is also the classical basis of the 
Tomonaga-Schwinger formulation of quantum field theory. 

Given any isolated system, containing any combination of particles, strings and fields 
and described by an action principle, one is lead to a reformulation of it as a parametrized 
Minkowski theory [6], with the extra bonus that the theory is already prepared for the cou- 
pling to gravity in its ADM formulation (Ref.[24]). This is done by coupling its action 
to an external gravitational field g^ v {x) and then by replacing the external 4- metric with 
Qab(j, <?)• In this way we get an action depending on the isolated system and on the em- 
bedding z^{t,(t) as configurational variables. Since the action is invariant under separate 
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r-reparametrizations and space-diffeomorphisms (frame-preserving diffeomorphisms [9]), ad- 
ditional first-class constraints are needed to ensure the independence of the description from 
the choice of the 3+1 splitting, namely from the convention chosen for the synchronization of 
distant clocks identifying the instantaneous 3-space to be used as Cauchy surface. The em- 
bedding configuration variables z^(t, a) are the gauge variables associated with this kind of 
special-relativistic general covariance and describe all the possible inertial effects compatible 
with special relativity. 

Let us come back to the discussion of free particles within the parametrized Minkowski 
theory approach. Since the intersection of a time-like world-line with a space-like hyper- 
surface, corresponding to a value r of the time parameter, is identified by 3 numbers a = fj{r) 
and not by four, each particle must have a well-defined sign of the energy. We cannot de- 
scribe, therefore, the two topologically disjoint branches of the mass hyperboloid simulta- 
neously as in the standard manifestly Lorentz-covariant approach. Then, there are no more 
mass-shell constraints. Each particle with a definite sign of the energy is described by the 
canonical coordinates ffi(r), Kj(r) (like in non-relativistic physics), while the 4-position of the 
particles is given by xf(r) = z^(t, ffi(r)). The 4-momenta pf(r) are ^-dependent solutions 
of pi — em\ = with the given sign of the energy. 

The system of N free scalar and positive energy particles is described by the action [6] 
S = J dTd 3 a£(T,a) = J drL(r), 

N 

C(t, a) = - 5 V ~ Vi{i~))mi\] 9tt(t, a) + 2g Tf (r, a)?)f(r) + g fS (r, d^fCr)^ (t), 



i=i 

N 



L ( r ) = -^misJg^T,^)) +2g Tf (r,fi i (T))fif(T)+gf S (T,ffi(T))^(T)fil(T), (3.4) 



i=l 



where the configuration variables are z M (r, a) and ffi(r), i=l,..,N. As said the action is 
invariant under frame-dependent diffeomorphisms [9]. In phase space the Dirac Hamiltonian 
is a linear combination of the following first-class constraints in strong involution 



T-(-^r,a) = p^(r,a) - l^(T,a)^25 3 (a - ffi(T))^ m? - 7^ (r, a) n if (r ) K is (r ) - 

i=i 

N 

Zf^r, a)Y s (r, a) <*V ~ ~ 0. (3.5) 

The conserved Poincare generators are 



p» = J d 3 ap^r, a), Jr = J d 3 a[z^r, a)p»(T, a) - z»(t, a)ff(r, a)}. 



(3.6) 
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B. The rest-frame instant form on the Wigner hyper-planes. 



Due to the special-relativistic general covariance implying the gauge equivalence of the 
descriptions associated to different admissible 3+1 splitting of Minkowski space-time in 
parametrized theories, the foliation can be restricted to space-like hyper-planes. In partic- 
ular, for each configuration of the isolated system with time-like 4-momentum, the leaves 
are best chosen as the hyper-planes orthogonal to the conserved total 4-momentum ( Wigner 
hyper-planes). Note that this special foliation is intrinsically determined by the configu- 
ration of the isolated system alone. This leads to the definition of the Wigner- covariant 
rest-frame instant form of dynamics [6], for every isolated system whose configurations have 
well-defined and finite Poincare generators with time- like total 4-momentum [8]. An inertial 
rest frame for the system is obtained by restricting the centroid world-line to a straight line 
orthogonal to Wigner hyper-planes. 

On a Wigner hyperplane, we obtain the following constraints [the only remnants of the 
constraints (3.5)] and the following Dirac Hamiltonian [6] 

N N 



e s - M sys « 0, k + = ^ Ki « 0, M sys = ^ y m i + «? » 

i=i i=i 

N 

H D = A(r)[e s - M sys ] - A(r) K h 



i=i 



±?(r) « -A(r)«"(p.) + e^u(p s ))X r (r), a£(r) = -A(t)«"(p.), 
<(r) = x » + u»(p s )T s + 6?(u(p s )) f dnX^). (3.7) 



The embedding describing Wigner hyper-planes is z^(r,a) = x%(t) + et^(u(p s ))a r , where 
e%(u(p)) = u^{p) = p^/^JTp* and ef{u{p)) = ( - u r (p); 5* - ! 4+i^r) are the columns of 
the standard Wigner boost, e%(u(p s )) = L^ u (p s ,p s ), connecting the time-like 4-vector p^ to 

its rest-frame form, p = s/ep 2 (1; 0). As a consequence the space indices "r" now transform 
as Wigner spin-1 3- vectors. 

— * 

The 3 Dirac's multipliers A(r) describe the classical zitterbewegung of the centroid x%(t): 
each gauge-fixing x( T ) ~ to the three first-class constraints k, + ~ gives a different 
determination of the multipliers A(r) and therefore identifies a different world-line x^P^{t) 
for the covariant non-canonical centroid. 

The various spin tensors and vectors evaluated with respect to the centroid are [6] 

Jg = XgP s — XgPg + Sg , 

Sr = [u»{p s y{u{p s )) - u\p s )e^u{p s ))]S: r + e^u(p s ))e v (u(p s ))S r s s , 
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N N 

Sf B = e A ^u{p s ))e B v {u{p s ))S^ = (s^ - ee - £ rf^m** + «*) , 

i=l i=l 
AT AT JV-1 

5 = 5 = X K, « ^j/ t X ^ - ?/ + X K + = ^p fl X 7f a . (3.8) 

i+1 i=l a=l 

Note that L^" = x^p v s —x u s p^ and S^" are not constants of the motion due to the classical 
zitterbewegung. 

Let us now summarize some relevant points of the rest-frame instant form on Wigner 
hyper-planes, since this formulation clarifies the role of the various relativistic centers of 
mass. This is a long standing problem which arose just after the foundation of special 
relativity in the first decade of the last century. It became soon clear that the problem of 
the relativistic center of mass is highly non-trivial: no definition can enjoy all the properties 
of the ordinary non-relativistic center of mass. See Refs.[38, 39, 40, 41, 42, 43] for a partial 
bibliography of all the existing attempts and Ref.[44] for reviews. 

Let us stress that, as said in the Introduction, our approach leads to a doubling of the 
usual concepts: there is an external viewpoint (see Subsection C) and an internal viewpoint 
(see Subsection D) with respect to the Wigner hyper-planes. While in Newton physics for 
any given absolute time there exists an absolute instantaneous 3-space and a unique center 
of mass with a unique associated 3- momentum P (vanishing in the rest frame), in special 
relativity the notion of instantaneous 3-space requires a convention for the synchronization 
of distant clocks. In other words we are forced: i) to introduce a 3+1 splitting of Minkowski 
space-time and an arbitrary accelerated observer; ii) for the sake of simplicity to restrict 
ourselves to inertial frames, namely to an inertial observer with a 4-momentum p^ and 
to foliations with space-like hyper-planes orthogonal to this 4-momentum (Wigner hyper- 
planes); iii) to describe the isolated system inside each of these instantaneous rest-frame 3- 
spaces, where its total 3-momentum vanishes. As a consequence, the unique non-relativistic 
center-of-mass 3-momentum P, vanishing in the rest frame, is replaced by two independent 

— * 

notions: i) an external 4-momentum p 1 ^ = (p°,p s = e s k s ) describing the orientation of the 
instantaneous 3-space with respect to an arbitrary reference inertial observer; ii) an internal 
3-momentum inside the instantaneous 3-space, which vanishes, k+ 0, as a definition of 
rest frame. The rest-frame instant form of dynamics is the natural framework to visualize 
this doubling. 

In the rest-frame instant form only four first-class constraints survive so that the original 
configurational variables z^ij, a), ffi(r) and their conjugate momenta p M (r, a), Kj(r) are 
reduced to: 

i) a decoupled point x%(t), p^ (the only remnant of the space-like hyper-surface) with 
a positive mass e s = \fep\ determined by the first-class constraint e s — M sys (M sys is 
the invariant mass of the isolated system). Its rest-frame Lorentz scalar time T s = is 
put equal to the mathematical time as a gauge fixing T s — r ss to the previous constraint. 
The unit time-like 4- vector u^{p s ) = p^/e s is orthogonal to the Wigner hyper-planes and 
describes their orientation in the chosen inertial frame. 

Here, x%(t) is a non-covariant canonical variable describing the decoupled canonical ex- 
ternal 4-center of mass. It plays the role of a kinematical external 4-center of mass and of 
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a decoupled observer with his parametrized clock (point particle clock). Its velocity x^(r) is 
parallel to p^, so that it has no classical zitterbewegung. The connection between the centroid 
x s ( r ) = z ^{ T i 0) an( i %s ( r ) an( i the associated decomposition of the angular momentum are 



x?(r) ^(t,£)=<(t)- 1 



e s (p° + e s ) 



« e 2 



V V/^s 1 V 

^ij = girfij'S?, Sf = s -^=, (3.9) 

P°s + V 6 P 2 s 

Now both Lf = x^p v s — x v s and Sf are conserved. 

ii) the 6N particle canonical variables ffi(r), /%(t) inside the Wigner hyper-planes. They 
are restricted by the three first-class constraints (the rest-frame conditions) k + = J2iLi ~ 
0. They are Wigner spin-1 3- vectors, like the coordinates a. 

The canonical variables x%, p^ for the external 4-center of mass, can be replaced by the 
canonical pairs [45] 



sr = 



rp Ps ' X s p s ■ X s I ^ -> /- Ps ~ g n 7* Ps / Q 1A \ 

-T s = = , e s = ±V e Ps> ^ = e s (x s -xj, fc s = — . (3.10) 

e s e s p° e s 

In the rest-frame instant form, this non-point canonical transformation can be summa- 
rized as 



x£ 




pt 







Z s 


Vi 


T 

-L S 


— * 

k s 


— * 



(3.11) 



The addition of the gauge-fixing T s — r m for the first-class constraint e s — M SJ/S w 0, 
building a pair of second-class constraints and implying A(r) = —1 (due to the explicit 
r-dependence of the gauge fixing), leads to the elimination of T s and e s . We are then left 
with a decoupled free point (point particle clock) of mass M sys and canonical 3-coordinates 
z s , k s . The position q s = z s /e s is the classical analogue of the Newton-Wigner 3-position 
operator [39] and shares with it the reduced covariance under the Euclidean subgroup of the 
Poincare group. 

The invariant mass M sys of the system, which is also its internal energy, replaces the 
non-relativistic Hamiltonian H rei for the relative degrees of freedom. For c — > oo we have 
M sys —mc 2 ^ H re [. M sys generates the evolution in the rest-frame Lorentz scalar time T s in 
the rest frame: in this way we get the same equations of motion as before the addition of the 
gauge fixing. This reminds of the frozen Hamilton- Jacobi theory, in which the time evolution 
can be reintroduced by using the energy generator of the Poincare group as Hamiltonian. As 
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a consequence, after the gauge fixing T s — r 0, the final Hamiltonian and the embedding 
of the Wigner hyperplane into Minkowski space-time become 



H D = M sys - A(r) • k + , 




= T 



(r) + e?(u(p s ))a r = a£(0) + U "(p 8 )r + e(f(u(p s )K, 




dr 



+ K(r), = ^(p s ) + e(f(«(p s ))A r (r), 



(3.12) 



where x^(0) is an arbitrary point. This equation visualizes the role of the Dirac multipliers 
as sources of the classical zittebewegung and consequently the gauge nature of this latter. Let 
us remark that the constant x£(0) [and 5^(0)] is arbitrary, reflecting the arbitrariness in the 
absolute location of the origin of the internal coordinates on each hyperplane in Minkowski 
space-time. 

Inside the Wigner hyperplane three degrees of freedom of the isolated system, describing 
an internal center-of-mass 3- variable <r com conjugate to k + (when the a com are canonical 
variables they are denoted q + ) are gauge variables. The natural gauge fixing in order to 
eliminate the three first-class constraints k + ~ is \ — 0+ ~ which implies A r (r) = 0: in 
this way the internal 3-center of mass gets located in the centroid x%(t) = z^ij.B = 0) of 
the Wigner hyperplane. 



C. The external Poincare' group and the external center-of-mass variables on a 
Wigner hyper-plane. 

The external viewpoint is proper to an arbitrary inertial Lorentz observer, describing the 
Wigner hyper-planes as leaves of a foliation of Minkowski space-time and determined by the 
time-like configurations of the isolated system. Therefore there is an external realization 
of the Poincare group, whose Lorentz transformations rotate the Wigner hyper-planes and 
induce a Wigner rotation of the 3- vectors inside each Wigner hyperplane. As a consequence, 
each such hyperplane inherits an induced internal Euclidean structure and an internal un- 
faithful Euclidean action , with an associated unfaithful internal realization of the Poincare' 
group (the internal translations are generated by the first-class constraints k + ~ 0, so that 
they are eliminable gauge variables), which will be described in the next Subsection. 

The external realization of the Poincare generators with non-fixed invariants ep 2 s = e 2 s ~ 

M; ys and W 2 = -ep 2 S s w -eM 2 ys S , is obtained from Eq.(3.8) (the four independent 
Hamiltonians [7] of this instant form, p° and J° l , are all functions only of the invariant mass 
M sys , which contains the possible mutual interactions among the particles): 



Js ~ X sPs X sPs ^ $s , 




M 2 ys +PZ = M, 



sys 




— * — * 

Ps = ^sk s ~ M S y S k, 
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N 



i=i 

e s + V e s + ^ ~t 

I Kir L.Sf-rsu Qu , fir„s_rsMC« 



(3.13) 



Given a canonical realization of the ten Poincare generators, one can build [18] three 
external 3- variables, the canonical 3-center of mass q s , the Moller 3-center of energy R s and 
the Fokker-Pryce 3-center of inertia Y s by using only these generators. For the rest-frame 
realization of the Poincare algebra given in Eqs.(3.13) we get 



R.i — 



P°s 



(X, 



^X°) 
P°s S) 



- _ ~ _ P±~o _ ££ _ D 

q s x s x s Jrig -f 

P°s e s 



S s X p s 

P°s(P°s + e s y 

S s xp s P°Rs + e s Y s 
P°s(P°s + e») Ps + ts 



Ys 



Qs + 



S s X ps 

e s (p°s + e s ) 



R s + 



S s X Ps 



Psts 



{R r s ,R s s } 



{P°s? " 

u:, q s s } = o, {y:,y:} 

— * — * — * 
Ps ■ Qs = Ps ' Rs = Ps " Yg = k s ■ z s , 



il s — J s R s x p S) 

\ c rsu gu _|_ S s ■ p s P s 

e s P° s L s ts(p°s + e s ) 

— » — * 

p s = => q s = Y s = R s . 



(3.14) 



All of these have the same velocity and coincide in the Lorentz rest frame where p s — e s (l; 0) 

Then, three external concepts of 4-center of mass can be defined (each having an internal 
3-location inside the Wigner hyper-planes) starting from the kinematics of the Wigner 
hyper-planes and from the above concepts of 3-centers of mass [38] : 

a) The external non-covariant canonical ^-center of mass x^ (with 3-location a), extension 
of the canonical 3-position vector q s (also named center of spin [42]). q s is the classical 
analogue of the Newton- Wigner position operator [39]. x 1 ^ is a frame- dependent pseu- 
dovector (q s does not satisfy the world line condition [38]), but it is canonical: {5^, x u s } = 0. 
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b) The external non-covariant and non-canonical M0ller J^-center of energy R% (with 
3-location <r#), extension of the M0ller 3-center of mass R s [40], which corresponds to 
the standard non-relativistic definition of center of mass of a system of particles with 
masses replaced by energies. R s does not satisfy the world line condition, so that is a 
frame-dependent pseudovector and moreover {R%, R v s } ^ 0. 

c) The external non-canonical but covariant Fokker-Pryce J^-center of inertia Y^ 1 (with 
3-location ay), extension of the Fokker-Pryce 3-center of inertia [41, 42]. Y^ 1 is a 4- 
vector by construction: it is the Lorentz transform of the rest-frame pseudo-world-line 
j^(rest)fi _ ^ est )^ Q f foe M0ller center of energy to an arbitrary frame. It holds {F/ 4 , Y"} ^ 0. 



Note that while the Fokker-Pryce Yf is the only 4- vector by construction, only x%(t) can 
be an adapted coordinate in a Hamiltonian treatment with Dirac constraints. 

To find the 3-locations on the Wigner hyper-planes, with respect to the centroid world- 
line x%(t) = z^ij, 0), of these quantities, we note [11] that, since we have T s = u(p s ) ■ x s = 
u(p s ) ■ x s = r on the Wigner hyperplane labeled by r, we can require Yf and R^ to have 
time components such that u{p s ) • Y s = u(p s ) ■ R s = T s = r. As a consequence, the 
following 3-locations are determined in Ref.fll]: a) for the world-line Y^{r) = z^(T,ay) = 

(x°(t); Y s (r)j we get a Y = R r + , with R + defined in Eq.(3.16) of next Subsection; b) for the 

pseudo-world-line x*{t) = z»(t,B) = (x°(t);x s (t)^ we get o r = a r Y + f+^f ; c) for the 

pseudo-world-line W a = z»(r, B R ) = (5 (r); i? s (r)) we get a r R = a Y + Hg^gg^l . 

It is seen that the external Fokker-Pryce non-canonical 4-center of inertia coincides with 
the centroid rrj,* + ^(r) carrying the internal 3-center of mass. 

In each Lorentz frame one has different pseudo- world-lines describing R% and x^\ the 
canonical 4-center of mass x 1 ^ lies in between Y^ and R 1 ^ in every (non rest)-frame. In an 
arbitrary Lorentz frame, the pseudo-world-lines associated with x^ and R^ fill a world-tube 
(see the book in Ref.[40]) around the world-line YJf of the covariant non-canonical Fokker- 
Pryce 4-center of inertia Y^. The invariant radius of the tube is p = V — eW 2 /p 2 = \S\/ y 1 ep 1 
where (W 2 = —ep 2 S 2 is the Pauli-Lubanski invariant when ep 2 > 0). This classical intrinsic 
radius delimitates the non-covariance effects (the pseudo-world-lines) of the canonical 4- 
center of mass x%. See Ref.[7] for a discussion of the properties of the Miller radius. At 
the quantum level p becomes the Compton wavelength of the isolated system times its spin 
eigenvalue y/ s(s + 1) , p i— > p = a/ s(s + Tjh/M = \/ s(s + 1)Xm with M = yfep 2 the 
invariant mass and Am — h/M its Compton wavelength. The critique of classical relativistic 
physics argued from quantum pair production concerns testing of distances where, due to 
the Lorentz signature of space-time, intrinsic classical covariance problems emerge: the 
canonical 4-center of mass x^ adapted to the first-class constraints of the system cannot be 
localized in a frame- independent way. Remember [6], finally, that p is also a remnant of 
the energy conditions of general relativity in flat Minkowski space-time: since the M0ller 
non-canonical, non-covariant 4-center of energy i? M has non-covariance properties localized 
inside the world-tube with radius p (see the book in [40]), it turns out that for an extended 
relativistic system with the material radius smaller of its intrinsic radius p one has: i) its 
peripheral rotation velocity can exceed the velocity of light; ii) its classical energy density 
cannot be positive definite everywhere in every frame. 
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D. The internal Poincare' group and the internal center- of- mass variables on a 
Wigner hyper-plane. 

Let us consider now the notions defined according to the internal viewpoint. They corre- 
spond to an unfaithful internal realization of the Poincare algebra: the internal 3-momentum 
k + vanishes due to the rest-frame conditions; the internal energy and angular momentum 
are given by the invariant mass M sys and by the external spin (angular momentum with 
respect to x^(t)) of the isolated system, respectively 

This internal realization of the Poincare algebra is built inside the Wigner hyperplane by 
using the expression of Sf B given by Eq.(3.8) (the invariants are M 2 ys — k 2 + pa M 2 ys > and 

W 2 = —e(M 2 ys — k\.)S s pa —eM 2 ys S s ; in the interacting case M sys and K are modified by 
the mutual interactions among the particles) 



N N 

M sys = H M = ^2\/m 2 + k 2 , K + = ^n i (« 0), 
i=i i=i 

N 1 

j = ffi x iu, r = s r = - € ruv s uv = s r s , 
i=i 

N i 

K = -^^m 2 + K 2 ff t , K r = J or = S7. (3.15) 
i=i 



As we shall see, K pa are the natural gauge fixings for the first-class constraints k + pa 0: 
this makes the internal realization even more unfaithful. 

In analogy with the external viewpoint, the determination of the three internal 3-center 
of mass can be achieved using again the group theoretical methods of Ref. [18]: 

i) a canonical internal 3-center of mass (or 3-center of spin) q + ; 

ii) a non-canonical internal M0ller 3-center of energy R + ; 

iii) a non-canonical internal Fokker-Pryce 3-center of inertia y + . 

Starting from the internal realization (3.15) of the Poincare algebra, we get the following 
Wigner spin 1 3-vectors: i) The internal Moller 3-center of energy R + and the associated 

— * 

spin vector Sr 



1 j? _ Ei=i V m2 i+z 2 m 



{Rl, <} = 6 rs , {Rl, M sys } = {i?;, R s + } = - j^e rsu S u R , 

m sys 1V1 sys 

{ST R , S s R } = e rsu (S u R - -^S R ■ k + <), {S r R , M sys } = 0. (3.16) 

1V1 sys 
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— * 

Note that the gauge fixing R + gives 



R + R + = {R + ,H D } = 



- A (r) ^ =1 « "Mr) « 0. (3.17) 



K 



k 



Furthermore, the internal boost generator of Eq.(3.15) may be rewritten as K = 
—M sys R + , so that R + ~ implies K « 0. 

ii) The canonical internal 3-center of mass q + and the associated spin vector S q [{q r + , q s + } = 
0, {q r + , k%} = 5 rs , {J r , q s + } = e rsu q u + , {S r q , S s q } = e rsu S^} 

^ J x ^ 

Q+ = r j 



^M% s -R\{M sys + JM% S - 4) 
K Jxk, 

+ 



+ 



— * 

K ■ K + K + 

M sys ^M^ ys -Kl[M sys + \J M^ ys - k\) ' 
« for k + « 0; {<f+, M SJ/S } = 



Sq — t/ — X — = + 

+ * x * + - ; /-« + « + ^ =j: (318) 



iii) The internal non-canonical Fokker-Pryce 3-center of inertia y + 



y + = g++ ; ^ =R + + 



£ | S q XK + _ MsyJh + sjM% s - & + y + 

M sys {M sys + ^M% s -Kl) M sys + ^M% s - 4 

{*;, y;> = — , 1 k + ? vj±^L l . 

M S ys^Mi ys - ^ L ^M^ s - 4(M sys + y/M^ - 4) J 



(3.19) 
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Note that on the Wigner hyper-planes, due to the rest-frame conditions k + ps 0, all the 
internal 3-centers of mass coincide, q + ~ R + ~ y + , and become essentially a unique gauge 
variable conjugate to k + . As a natural gauge fixing for the rest-frame conditions k + ps 0, 
we can add the vanishing of the internal Lorentz boosts K = —M sys R + ps 0: this implies 
A(r) ~ and is equivalent to locate the internal canonical 3-center of mass q + in a = 0, i.e. 
in the external centroid x^(t) = z^(t, 0) origin of the internal 3-coordinates in each Wigner 
hyperplane. Remember that the centroid x%(t) corresponds to the unique special- relativistic 
center-of-mass-like world-line of Refs.[46]. With these gauge fixings and with T s — r ps 0, the 
world-line x^(r) of the centroid becomes uniquely determined except for the arbitrariness in 
the choice of x^(0) 

X f+)^ T = Ts ) = a£(0) + u»(p s )T s , (3.20) 

and coincides with the external covariant non-canonical Fokker-Pryce 4-center of inertia, 
x%(t) = x^(0) + Yf. It can also be shown that the centroid x^\r) coincides with the Dixon 
center of mass of an extended object [29] as well as with the Pirani [47] and the Tulczyjew 
[48] centroids. 

— * 

Note that in the non- relativistic limit all the quantities q + , R + , y + tends the the non- 
relativistic center of mass q nr = m%r>% . 

We are left with the problem of the construction of a canonical transformation bringing 

from the basis ffi, K, iy to a new canonical basis q + , 0), p q>a , n q:a , in which S q = 

\ - v i . 

rfi_ 

Let us stress that this cannot be a point transformation, because of the momentum de- 
pendence of the relativistic internal 3-center of mass q + . The canonical transformation (3.21) 
will be constructed in the next Subsection by using the method of Gartenhaus-Schwartz [49] 
as delineated in Ref.[50]. 

In conclusion, at the relativistic level the non-relativistic Abelian translation symmetry 
generating the non-relativistic Noether constants P = const, gets split into the two follow- 
ing symmetries: i) the external Abelian translation symmetry whose Noether constants of 
motion are p s = e s k s ~ M sys k s = const, (its conjugate variable being the external 3-center 
of mass z s ); ii) the internal Abelian gauge symmetry generating the three first-class con- 
straints k + ps (rest-frame conditions) inside the Wigner hyperplane (the conjugate gauge 
variable being the internal 3-center of mass q + ps R + ps y + ). Of course, its non-relativistic 

— * 

counterpart is the non-relativistic rest-frame condition P ps 0. 

E. The canonical transformation to the internal center-of-mass and relative vari- 
ables for N free particles. 

Since q + and k + are known, we have only to find the internal conjugate variables appearing 
in the canonical transformation (3.21). They have been determined in Ref. [11] by using 
the technique of Ref. [49] and starting from a set of canonical variables defined in Ref. [6]. 





Pqa 


«+ 


Kqa 



(3.21) 



29 



Precisely, starting from the naive internal center-of-mass variable ff + = Y2iLi Vh we applied 
with definition of relative variables p a , ir a based on the following family of point canonical 
transformations [the numerical parameters 7 ai satisfy the relations in Eqs.(2.2)] 



v+ 


Pa 







a = 1,..,N- 1, 



V+ 



Pa = 



K; 



I ""I 

v a=l 
^ AT N 

i=l i=l 

AT AT 



= —K+ + ^^7™^, 



o=l 



i=l v i=l 



{Pa, <} = Sa b S T 



(3.22) 



Then, (in Appendix B of Ref. [12]), we gave the closed form of the canonical transformation 
for arbitrary N, which turned out to be point in the momenta but, unlike the non-relativistic 
case, non-point in the configurational variables. 

Explicitly, for N = 2 we have 



M, 



sys 



7T„ 



Pa 



^m\ + k\ ffi + y/ml + k- 



•2 -» 
2^2 



(?7l X Ki + 772 X K 2 ) X «+ 



-*2 
«1 



+ 



- 4 (M SJ/S + y/M^ - 4) 



(V^l + «1 m + V m 2 + ^2 %) • «+ ^ 



/?! + K 2 ~ 0, 
7T — 



M% s - 4 (M SJ/S + JM% S - k 



m\ + k\ — J ml + k\) — 



+ 



K + ■ IT 



(M sys - X /M? ys - k 



7T, 



y/ml + Kl sJm\ + K 2 



i) 



■y/m^ + ifl y/mf + vff M sys ^M% s -k\ 
M sys = sJm 2 + k\ &M = sjm\ + i?l + sjm 2 2 + itl. 



(3.23) 



The inverse canonical transformation is 
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Vi = Q+ 



^M 2 + & + (M + ^M 2 + k\ 



1 

+ 2 



2 M 7i q ■ k + + {m\ - m 2 ,) a/-^ 2 + k 
M 2 ^/M 2 + 



Pq 



Pq ■ K + Kg 



M y/M 2 + K 2 + 



<?+ + 



1 r 



vi+i m i - m l 



.M 2 



7T, •«+[!- 



A* 



-2 .M y/Af 2 - //: v ' 
+ K " m\) Jm 2 + kI)] k + + (-)* +1 



, . v M 2 + 4 - A4)] + 



i+1 7T 



(3.24) 



For N > 2, the Hamiltonian M sys = Xlili \/ m l + N Z)^ -1 lailu^qa • n qb for the 
relative motions in the rest frame instant form, is a sum of square roots each containing a 
(AT — 1) x (N — 1) matrix K7} b = ^7^7^ = Kzh [note that in the non-relativistic limit 

only one such matrix appears, namely = YliLi ^ K-uLb °^ Eqs.(2.3)]. The existence 
of relativistic normal Jacobi coordinates would require the simultaneous diagonalization of 
these N matrices. This is impossible, however, because 

[^(i) 1 ' K (j)1 ab = G(ij)ab = ~G{ij)ba = -G(ji) ab = —N[j ai J bj - Jajlbi]- (3.25) 

There are ^N(N — 1) matrices Gij, each one with j(N — 1)(N — 2) independent elements. 
While the conditions G^ ab = are jN(N — 1) 2 (N — 2), the free parameters at our disposal 
in the are only | (AT — 1) (AT — 2). For N=3, there are 3 conditions and only 1 parameter; 
for N=4, 18 conditions and 3 parameters. 

In conclusion: it is impossible to diagonalize the N quadratic forms under the square 
roots simultaneously, no relativistic normal Jacobi coordinates exist, and reduced masses 
and inertia tensor cannot be defined. 

The relative Lagrangian can be worked out in the special case of N=2 with equal masses 
(mi = m 2 — m). It results 

LrMfi = -m^A-j?. (3.26) 

so that the relative velocity is bounded by |p| < 2. 

Let us write p = pp with p — \p\ and p = With a single relative variable the 
three Euler angles 8 a are redundant: there are only two independent angles, identifying the 
position of the unit 3- vector p on S 2 . We shall use the parametrization (Euler angles 6 1 = 0, 

02 = 0,03 = 0) 
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p r = R rs {9A)K = (RMRyi&y' 'ft, Po = (0,0,1). (3.27) 

In analogy with Subsection IA of the non-relativistic case, we get the following body 
frame velocity and angular velocity (p is the only shape variable for N=2) 



& = R T ™p° = p{R T RY s pl + ppl = pe"*u>" + ppl = p(u x p o y + p%, 

- = A e urs (RTxysj = ^1 = _ sin ^2 = ^ Q), 

v 2 = I(p)u 2 + p 2 , I(p)=p 2 . (3.28) 

The non-relativistic inertia tensor of the dipole I nr = pp 2 ( p = ^^m, * s ^ ne re duced 
mass) is replaced by / = I nr /p = p 2 ■ The Lagrangian in anholonomic variables become 



L{u, p, p) = - m ^4-I(p)uj 2 -p 2 . (3.29) 

It is clear that the bound \p\ < 2 puts upper limitations upon the kinetic energy of both the 
rotational and vibrational motions. 
The canonical momenta are 

- dL ml(p)u) dL mp 

b = — = , 7T = — = - (3.30) 

0UJ <U-I{p)uJ 2 -p 2 ° P J 4 - I{p)u 2 - p 2 



Note that there is no more a linear relation between spin and angular velocity. 

When |p| varies between and 2 the momenta vary between and oo, so that in phase 
space there is no bound from the limiting light velocity. This shows once more that in 
special relativity it is convenient to work in the Hamiltonian framework avoiding relative 
and angular velocities. 

Since we have \/4 — I(p)uj 2 — p 2 = , 2m = , the inversion formulas become 

V V ' ' r ^m>+I-^{p)S 2 +-K- 2 

2I-\p)S 



uj = - I{p)u 2 -p 2 = — , 

^m 2 + I^(p)S 2 + n 2 

p = y4-I{p)Z 2 -p 2 = 2n = =. (3.31) 

m ^m 2 + I- 1 (p)S 2 + 7i 2 



Then, the Hamiltonian results 



M sys = rrp + S-u-L = 2\Jm 2 + I-^pjS 2 + vr 2 . (3.32) 
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F. Dynamical body frames and canonical spin bases for N relativistic particles. 



For isolated systems the constraint manifold [7] is a stratified manifold with each stra- 
tum corresponding to a type of Poincare orbit. The main stratum (dense in the constraint 
manifold) corresponds to all configurations of the isolated system belonging to time-like 
Poincare orbits with ep 2 ~ eM^ > 0. As said in Ref. [20], this implies that the inter- 
nal 3-center-of-mass coordinates have been adapted to the co-adjoint orbits of the internal 
realization Poincare group. But, since the second Poincare invariant (the Pauli-Lubanski 
invariant W 2 = —p 2 S 2 s ) does not appear among the canonical variables, this canonical basis 
is not adapted, as yet, to a typical form of canonical action of the Poincare group [18] on the 

phase space of the isolated system. However, remember that the scheme A for the internal 

s 2 

realization of the Poincare group [18] contains the canonical pairs k + , q+, S^, arctg ^r, and 
the two Casimirs invariants \S q \ — ^J—W 2 /M^ ys , M sys = ^2f =1 y/m 2 + k 2 . 

As a consequence, as shown in Ref. [20], it is possible to construct a canonical basis 
including both Poincare invariants in such a way that all of the internal coordinates are 
adapted to the co-adjoint action of the group and the new internal relative variables are 
thereby adapted to the SO (3) group. 

In conclusion, in the rest-frame instant form of dynamics the construction of the internal 
dynamical body frames and of the internal canonical spin bases is identical to that of the 
non-relativistic case. Only the form of the relative Hamiltonian, the invariant mass M sys , is 
modified. 



G. Action-at-a-distance interacting particles and relativistic orbit theory. 

As shown in Ref. [6] and its bibliography (see also Ref. [51]), the action-at-a-distance in- 
teractions inside the Wigner hyperplane may be introduced under either the square roots 
(scalar and vector potentials) or outside (scalar potential like the Coulomb one) appearing 
in the free Hamiltonian. Since a Lagrangian density in presence of action-at-a-distance mu- 
tual interactions is not known and since we are working in an instant form of dynamics, the 
potentials in the constraints restricted to hyper-planes must be introduced by hand. The 
only restriction is that the Poisson brackets of the modified constraints must generate the 
same algebra of the free ones. 

In the rest-frame instant form the most general Hamiltonian with action-at-a-distance 
interactions is 



M syS;int = Y t y/™? + U i + [h -W + V, (3.33) 

i=i 

where U = U(K k ,ff h - ff k ), V { = V^kj^fji - rfj#), V = V (\ffi - ffj\) + V ' {k h fji - ffj). 

If we use the canonical transformation (3.21) defining the relativistic canonical internal 3- 
center of mass (now it is interaction-dependent, q^ 1 ^) and relative variables on the Wigner 
hyperplane, with the rest-frame conditions k + k 0, the rest frame Hamiltonian for the 
relative motion becomes 
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JV-l 

m* + U i + [VNj2larf qa - Vtf + V, (3.34) 

a=l 

where 

JV-l ^ JV-l 

Ui = U([VN^-f ak 7r qa ,-=y2(-f ah --f ak )p qa ), 
— 7 V -/V — i 

a=l v a=l 

JV-l JV-l 

a=l * a=l 

^ JV-l JV-l JV-l 

^ = Vo (l7^E( 7a *~ 7a ^ a l) + v '([^X^^^ 

v a=l a=l * a=l 

In order to build a realization of the internal Poincare' group, besides M sys ^ nt we need to 
know the potentials appearing in the boosts K~i nt (being an instant form, k + pa and J are 
the free ones). We need therefore the rest-frame energy- momentum tensor of the isolated 
system (see later). 

Since the 3-centers R + and q + become interaction dependent, the final canonical basis q + , 
k+, p q ai TTqa is not explicitly known the interacting case. For an isolated system, however, 
we have M sys = a/-M 2 + k 2 + pa M with M independent of q + ({M sys , k + } = in the 
internal Poincare' algebra). This suggests that the same result should hold true even in the 
interacting case. Indeed, by its definition, the Gartenhaus-Schwartz transformation gives 
Pqa ~ Pa, ^qa ~ also in presence of interactions, so that we get 



JV 



i=l \ 



M sysMt \ R+=0 = (^ yjm* + U t + (Ki - Vtf + v) U +=0 = y] Ml nt + = 

i 

= M int \z + =o = \Jmt + U, t + (K t -V^ + V, (3.36) 

i 

where the potentials Ui, V ^ V are now functions of ir qa • n qb , 7r qa ■ p qb , p qa ■ p qb . 

Unlike in the non-relativistic case, the canonical transformation (3.23) is now interaction 
dependent (and no more point in the momenta), since q + is determined by a set of Poincare' 
generators depending on the interactions. The only thing to do in the generic situation is 
therefore to use the free relative variables (3.23) even in the interacting case. We cannot 
impose anymore, however, the natural gauge fixings q + ~ (K pa 0) of the free case, since it 
is replaced by g+ n ^ ~ (namely by Ki n t ~ 0), the only gauge fixing identifying the centroid 
with the external Fokker-Pryce 4-center of inertia also in the interacting case. Once written 
in terms of the canonical variables (3.23) of the free case, these equations can be solved 
for q + , which takes a form q + pa f(p aq , ir aq ) as a consequence of the potentials appearing in 
the boosts. Therefore, for N = 2, the reconstruction of the relativistic orbit by means of 
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Eqs.(3.24) in terms of the relative motion is given by (similar equations hold for arbitrary 
N) 



fji(r) « q+ipq^q) + - 

<W = 4-(^.«)=^(0) + ^(P^ + <(^))^W' (3-37) 

While the potentials in M sys ^ int determine p q (r) and 7r 9 (r) through Hamilton equations, 
the potentials in K int determine q+(p q , n q ). It is seen, therefore - as it should be expected - 
that the relativistic theory of orbits is much more complicated than in the non-relativistic 
case, where the absolute orbits ffi(t) are proportional to the relative orbit p q {t). 

A relevant example of this type of isolated system has been studied in the second paper 
of Ref.[6] starting from the isolated system of iV charged positive-energy particles (with 
Grassmann- valued electric charges Qi = 9* 9, Q\ = 0, QiQj = Qj Qi ^ for i ^ j) plus 
the electro-magnetic field. After a Shanmugadhasan canonical transformation, this system 
can be expressed only in terms of transverse Dirac observables corresponding to a radiation 
gauge for the electro-magnetic field. The expression of the energy-momentum tensor in this 
gauge will be shown in the next Subsection IVC (where K int can be calculated). In the 
semi-classical approximation of the second paper of Ref.[6], the electro-magnetic degrees 
of freedom are re-expressed in terms of the particle variables by means of the Lienard- 
Wiechert solution, in the framework of the rest-frame instant form. In this way the exact 
semi-classical relativistic form of the action-at-a-distance Darwin potential in the reduced 
phase space of the particles has been obtained. Note that this form is independent of the 
choice of the Green function in the Lienard-Wiechert solution. In the second paper of Ref. [6] 
the associated energy-momentum tensor for the case N = 2 [Eqs.(6.48)] is also given. The 
internal energy is M sys = ^ M 2 + k\ M = Ya=i V m i + 7?2 + ^ 2 [i + V^TT 2 , 7?-^)] where 
V is given in Eqs.(6.34), (6.35) [in Eqs. (6.36), (6.37) for mi = 777.2]. The internal boost 
K int [Eq.(6.46)] allows the determination of the 3-center of energy R + = —jf^ ~ <f+ ~ y+ 
in the present interacting case. 



Pq- 



■ +1 mi - m 2 
/ ™ Pi -1 



m, 



H. An Example of Deformable Continuous System: the Classical Klein-Gordon 
Field. 

Consider now the problem of separating the relativistic center of mass for isolated special- 
relativistic extended systems. In Ref. [46], mainly devoted to the same problem in general 
relativity, it is shown that in special relativity there is a unique world-line describing this 
notion, corresponding to the centroid x%(t) of the rest-frame instant form. The first attempt 
to define a collective variable for a relativistic extended body in the Hamiltonian formulation 
was done in Ref. [25] in the case of a configuration of the Klein-Gordon field (which is used, 
e.g., in the description of bosonic stars). Let us show how the previous kinematical formalism 
is working in this case. 

The rest-frame instant form for a classical real free Klein-Gordon field [26] 0(r, a) = 
0(z m (t, a)) (4>(x) is the standard field) on the Wigner hyper-planes is built starting from 
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the reformulation of Klein-Gordon action as a parametrized Minkowski theory. The relevant 
quantities, as well as the internal and external Poincare' generators, can be worked out as in 
Section III, Subsection A, B, C and D, for N free relativistic particles [q A = (q T = u(q); q r ), 



= y/m 2 + q 2 , Q{q) = (2vr) 3 2u{q), dq = d 3 q/Q{q)] 



(r, a) 
7T (r, a) 



a(r,q) — / d a [u (q) (j) (r, a) + in (r, a)] e 
a* (r, q) 



i(u)(q)r-q-a) 



J d 3 a [u (q) (r, a) - in (r, a)} e^M^"**), 

J dq [a (r, q) e ~i(^-^) + a * (r, q) e +*Mfr-**)] , 
= -i J dqu (q) [a (r, q) e^M^-**) _ a * ( r> e +*(^>-<f-a)] ^ 

{<f>(ra), tt(t, 5i} = 5 3 (<? - a ± ), [a (r, g) , a* (r, fc) } = -in (q) 5 3 (q- k) . 

(3.38) 



Ma, 



y dV vr 2 + (<90) 2 + m 2 2 (r, a) = J dq u (q) a* (r, q)a(r,q), 



d 3 a[ird(f)](T, a) = / dqqa* (r : q) a {r,q) ~ 0, 



-i J dqa* (r, g) (q r - q s -^j a (r, q) , 

-rP; + 1 y dV a' [vr 2 + («90) 2 + m 2 2 ] (r, a) 

f ~ d 
~ rP l + l J dqu(q)a* (r,q) —a(r,q) , 



ss = j; s |^ =0 = y d 3 a{a r [nd s <t>](r,a) - (r <-> s)}\ p ^ =Q . 



(3.39) 



Here is the usual canonical non-covariant 4-center of mass of Eq.(3.9). We are working 
on Wigner hyper-planes with r = T s , so that the Hamiltonian (3.12) is H D = M^ — A(r) - Pa,. 
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We want to construct four variables X£[<f>, it] = (X^X^) canonically conjugated to 

P^[0,7r] = (P^;Pej,). First of all we make a canonical transformation to modulus-phase 
canonical variables 



a*(r,q) = y/TJ^je^^, 



I(r,q) = a* (r, q)a(r,q), 



V (r, q) = — In 



a (t, q) 
a* (r, q) 



{I(r,q)^(r,q')} = n( q )5 3 (q-q*) 
In terms of the original canonical variables 0, ir, we have 



(3.40) 



I(r,q) = Jd 3 aJ rfV Kg)0(r, a) - m(r, a)} [u(q)<j>(T, a') + tn(r,a% 

, 1, r fd 3 a[uj(qU(T,a)+in(r,a)}e i ^ T -^ 

= 7T ln 



/ dV [u;(g)0(r, <?') - i7r(r, B')]e- i ^ T -^ s ') 



/ d 3 a [w(g)0(r, a) + Z7r(r, a)] e ~^" 



/ dV [w(g)0(r, o") - i7r(r, ^Je^' J ' 



(r, a) = j dq yj I (r, g) [ e M^-) -^tor-fa) + e -*v(r,S)-H(a;( 9 )r-^j ^ 
7T ( r , «?) = -i y dgcu (g) v^MO [ e M^-)-*(u;(?)T-<f-a) _ e -<iP(r,«)+i( w (?)r-f.9)] _ ( 3 41 ) 
Then, the internal Poincare charges of the field configuration take the form 



dqu(q)I (r,q) , P = / dqql (r, g) w 0, 
dg/ (r, g) (g r ^ - g s A ) ^ (r, g) , 



jrs _ 

* ~ J ' 

J T i = ~rP;- Jdqcv(q) 
while the classical analogue of the occupation number is [A = — d 2 } 



d 

I (r, q) g-^V (r, q) , 



N$ = J dqa*(r,q)a(r,q) = J dql(r,q) = 

= \f d 3 a[n 1 ■ + v /m 2 + A]0(r,o i ). 
* J ^/m z + A 



(3.42) 



(3.43) 
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1. Definition of the collective variables. 
Define the four functional of the phases 



X; = J dqu(q)F T (q) up (r, q) , = J dqqF (q) up (r, q) , 

^{x;,x;} = o, {x;,x;} = o. (3.44) 

depending on two Lorentz scalar functions F T (q), F(q). Their form will be restricted by the 
following requirements implying that X^ and P^ are canonical variables 

{p;,x;} = 1, {p;,x|} = -5", {p;,x;} = o, {p;,x;} = o, (3.45) 

Since {P^,X^} = j dqu 2 (q) F T (q) and {P^,X|} = j dqq r q s F(q), we must impose the 
following normalizations for F T (q), F(q) 



J dqcu 2 (q)F T (q) = 1, J dqq r q s F (q) = -5 T 



(3.46) 



Moreover, {P£,X;} = Jdquj(q)q r F T (q) and {P^X;} = / dquo (q) f F (q) , imply the 
conditions 

J dqu(q) q r F T (q) = 0, J dqu(q) q r F (q) = 0. (3.47) 

that are automatically satisfied since F T (q), F(q), q = \q\, are even under q r — > — q r . 
A solution of Eqs.(3.47) is 

™-/ x 167T 2 4tt 2 . 487T 2 / Att „2 , . 

r ^ rrTrf^ ' F(g) = 4^ 2 + 9 2 e-^*. (3.48) 

mq z ^m z + q z rnq 



The singularity in g = requires </?(r, q) ^ q ^o q v , rj > for the existence of X^, X^. 

Note that for field configurations 0(r, <r) such that the Fourier transform 0(r, q) has 
compact support in a sphere centered at <f = of volume V, we get X^ = J ^^^(t, q), 

X^^JdVMr,q). 



2. Auxiliary relative variables. 



As in Ref. [25], let us define an auxiliary relative action variable and an auxiliary relative 
phase variable 



I (r, q) = I (r, q) - F T (q) P^cu (q) + F (q) q- P , 

(r, q) = <p (r, q) - to (q) X% + q- X*. (3.49) 
The previous canonicity conditions on F T (q), F(q), imply 
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dqu (g) I (r, q) =0, J dqq* I (r, q) = 0, 

J dqF T (q) co (g) (r, q) dq = 0, | dqF (q) q l (r, q) = 0. (3.50) 
Such auxiliary variables have the following non-zero Poisson brackets 

{/(t,£),0(t,$)} = a(£,$) = 

= (A;) 5 3 (k-q^j- F T (k) u (k) u (q) + F (k) k ■ q. (3.51) 

The distribution A (k,q) has the semigroup property and satisfies the four constraints 



dqA (k, g) A (q, k^j = A (jfe, fc') , 
J dqu (g) A (q, fc) = 0, | dgg r A (q, k) = 0, 

y dqF T (q)uj(q)A(k,q) = 0, J dqq r F (q) A (k,qj = 0. (3.52) 

At this stage the canonical variables J(t, g), </?(r, g) for the Klein-Gordon field are replaced 
by the non-canonical set , P^ , X^ P^, I(r,q), 0{r,q) with Poisson brackets 

{p;,x;} = i, {p;,x;} = -5 rs , {p;,x;} = o, {p;,x;} = o, 
{x;,x;} = o, {x;,x;} = o, {pf,p*} = o, a,b = (t,t), 

{p;, / (r, g)} = 0, {p£, / (r, g)} = 0, { J (r, g) , X;} = 0, { J (r, g) , X^} = 0, 
{X;,0(r,q)} = O, {X;,0(r,q)}=O, {P£, (r, g)} = 0, {P^, (r, g)} = 0, 
{/(r,jfe),^(r,g)} = ft (A;) 5 3 - g) - P r (Jfe) (Jfe) (g) + F (k) k ■ q. (3.53) 

Note finally that the generators of the internal Lorentz group are already decomposed 
into the collective and the relative parts, each satisfying the Lorentz algebra and having 
vanishing mutual Poisson brackets 

jrs r rs i qrs 

• J 4> ^ (p ' 

^ S = ^ - X; P ;, S? = J dqi(r, q) A - q^ (r, q) , 

jrr t tv I qtt' 

^7 = - r\P; - X;P;, S; r = -Jdqu (g) J (r, g) ( T) g1 . (3.54) 
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3. Canonical relative variables. 



Now we must find the relative canonical variables hidden inside the auxiliary ones. They 
are not free but satisfy Eqs.(3.50). As in Ref.[25], let us introduce the following differential 
operator [A LB is the Laplace-Beltrami operator of the mass shell sub-manifold H\ (see 
Ref.[25, 26])] 



£>„- 



3 - m A LB 

3 



o 2 

= 6 — m 



" ( d "\ ^ 2 c) 1 / ° d 

^ \ dq i J m 2 ^ ^ So* m 2 \ ^ ^ <9g* 
i=i x ' i=i \ «=i 



(3.55) 



Note that, being invariant under Wigner's rotations, is a scalar on the Wigner hyperplane. 
Since u (q) and q are null modes of this operator [25], we can put 



I(T,q)=VfH(r,q), H(r,q) — j dkQ (q,k) I (r,k) . 
with Q (^q, itj the Green function of (see Refs. [25] for its expression) 



T>^3{q,k) =n(k)5 3 (fc-g). 



Like in Ref.[25], for each zero mode f (q) of V$ \V^f (q) 
J dq f (q)I( T ^\ < °°> integrating by parts, we get 



(3.56) 

(3.57) 
0] for which 



J dq f (q)i(r,q) = J dq f (q)V ? U(r, g) 
1 f o 9 /m 2 5 rs + g r g 



2(2tt) 3 J * <9g r V 
The boundary conditions (ensuring finite Poincare generators) 



/o(5)^H(r,5)-H(r,5)A /o(5) 



(3.58) 



H(r, g) 
H(r, g) 



+ o g 



-l+e 

J 

-3-ct 



e > 0, 
a > 0, 



-q^oo 9 " o- > u, (3.59) 

imply f dq f (q)I(T,q) = (so that the first two conditions (3.50) are also satisfied), or 

dqf (q)I(r,q) = P; J dq u (q) f (q) F T \q) - P • J dqqf (q)F(q). (3.60) 

It is shown in Ref.[25] that, by restricting ourselves to field configurations for which 
I{j,q) ^ q ^o(l 3+ri with i] G (0, 1] and by imposing the following restriction on <p(r,a) and 
7t(t, a) = d T (f)(r, a) 
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Pio = const. I dqq l 2 F^— ^; -q 2 )Y l0 (9, <p) 



3 
2 ! 



d 3 a / dV'e^P- ) 



K + g^^^^^^a^+Tr^a)^,^) =0. (3.61) 



we can identify the class of the Klein-Gordon field configurations that is compatible with 
the previous canonical transformation and lead to a unique realization of the Poincare group 
without any ambiguity 

We can satisfy the constraints (3.50) on <p(r, q) with the definition [D$u(q) = V^q = 0] 



(p(r,q) = J dk J dk'K(r,k)g(k,k')A(k',q), 
K(r,q) = V<fp(T,q) =V ( fp(T,q), 



^oog 1 e , e > 0, 



which also imply 



{H (r, q) , X;} = 0, {H (r, q) , P;} = 0, 

{H(T,q),X;} = 0, {H(r,q),P;}=0, 

{K (r, q) , X;} = 0, {K (r, q) , P;} = 0, 

{K(r,g),X;}=0, {K (r, q) , P£} — 0, 

{H (r, g) , K (r, ?)} = (g) 5 3 (q- <?) . 

The final decomposition of the internal Lorentz generators is 



(3.62) 



(3.63) 



jrs r rs , Qrs 

Lrs -v-r r>s r>r 

<t> — A 0-W» ~~ ^V4> 

JTT T TT I orr 



9 



<9A; S 



K 



(r, fc) , 



= (x; - r)p; - x;p;, s; r = -J dqu ( q ) h (r, $ A K ( r , $) . (3. 



64) 



4- Field variables in terms of collective-relative variables. 
We have found the canonical transformation 



I(r,q) = F T (q)u(q)P; - F(q)q- P^ + V <r H(r,q), 

<p(r, q) = J dk J dk'K(r, k)G (X fc') A g) + cu (q) X\ - q- 
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= P; j dqu(q)F^q) - ■ J dqqF(q) + J dqV^H(r, q) = 



j dqV ? U(T,q), 



m 
c = c 



(m) =m J dqu(q)F T (q) = 2^ 



;/t/ = 2e 4 - 



\Jm 2 + q 2 

with the two functions F T (q) ) F(q) given in Eqs.(3.48). Its inverse is 



\Jx 2 — 



-e m 



(3.65) 



PI = J dqcu(q)I(r, q)= l -J d 3 a [n 2 + (<90) 2 + m 2 2 ] (r, a), 

P<t> = J dqql(r,q) = J d 3 a 7r<90 (r, <?), 

X; = J dqco(q)F T (qMr,q) = r + 

u(q) J d 3 ae^ ff (f)(T, a)+ij dW^T^r, a) 



+ 



def 



4.7V 2 



d 3 q- 



-In 



2mm J q 2 ^(q) >-w(q) J d 3 ae l ^ ff (j)(T,a) —if d 3 ae 1 & 3, k{t,&) 



r + XL 



a — JV±r±. — J\.±r± 



X<t> = J dqqF(q)(p(T,q) = 



- d 3 q^-e ^ q In 
nj V 



^m 2 + q 2 J d 3 ae iq - ff (f)(T, a) + i f dW^^T, a) 
^m 2 + q 2 j d^ae-^cpir, a) - i J d 3 ae~ i ^n(r, a) 



H(r, q) = j dkg(q,k)[I(r,k)-F T (k)uj(k) / ^(^/(r, &) + 
+ F(k)k- J dq&Ifafr)] = 



d 3 o 1 d 3 o 2 |tt(t, <7i)tt(t, a 2 ) J dkQ(q, k) J dkiA(k,k 1 )e ikl '^ 1 ~ ff2 ^ + 
+ 4>{r,a 1 )4>{r,a 2 ) J dkG(q,k) J dk 1 u 2 (k 1 )A(k,k 1 )e ikl ' {3l ~ S2) - 
- i (n(r, CTi)0(r, a 2 ) + 7t(t, <t 2 )0(t, 

dkG(q,k) [ dhujik^Aik.k^e^ 1 -^], 



K(r,q) = V<?p(T, q) = V#p(t, q) 



/ d 3 a [u (q) (r, a) + ivr (r, a)] e _< ** 
JT dV [u (q) (r, a') - in (r, a')] e*^ 



(3.66) 
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The remaining canonical variables a(r,q), <f>(r,a), 7t(t, a), can be worked out in terms of 
the final ones 



a(r,q) = ^F-(q)u(q)P; - F(q)q- P^ + P f H(r, q) 

e i[u>(q)X^-q-X^]+iJdkJdk'K(T,k)g(k,k')A(k',q) 



pr 



0(r, a) 



J dq^FT{q)u{q)Pl - F(q)q- P + %H(r, q) 

e -i[ W ( g )(r-Xj)-g-(a-X )]+i/dfc/dfe'K(r,fc)6(fe,fc')A(fe',q) + 
+ e i[ W ( g )(r-X;)-g-(a ! -X )]-i/dfe/ ( ife'K(r,fc)g(fc,fc')A(fc',g)j = 

2 / dgA^r; Pf , H] cos \q- a + B^r; X* K] 



(3.67) 



7r(r,5) = -ij dqco(q)y/F-(q)u;(q)P; - F(q)q- P^ + P ? H(r, g) 

e -i[ W (?)(r-Xj)-g-(o ! -X^)]+i/ C ifc/dfe'K(T,fe)£;(fe,fe')A(fe',q) _ 
_ e +i[o;( 9 )(r-Xj)-q.(5 ! -X^)]-i/dfe/dfe'K(T,fc)£;(fc,fc')A(fe',q) 

= -2^ dgw(g)A^T;P/,H]sin [<f • a + B ? -(r; X$, K] 

A q ir; Pf, H] = ^F^q)uj(q)P; - F(q)q- P^ + ^H(r, q) = y/l{ fjj, 
B^r-X^K] = -q.X (f) -Lu(q)(r-X T 4> ) + J dkdk'K(T,k)g(k,k')A(k',q} = 

= <p(r,q) - w(q)r. (3.68) 
Summarizing, the Klein-Gordon field configuration is described by: 

i) its energy P£ = M^, i.e. its invariant mass, and the conjugate field time-variable 

def ~ ~ 

X^ = t + X^ (r = T s ), which is equal to r plus some kind of internal time X^ (note that in 
the N-body case this variable does not exist: M sys is a given function of the other canonical 
variables and not an independent canonical variable like here); 

ii) the conjugate reduced canonical variables of a free point X^, P^ ps (see fj + and k + 
of Eqs.(3.22) in the N-body case); 

iii) an infinite set of canonically conjugate relative variables H(r, q), K(r, q) (p a and 7? a of 
Eq.(3.22) in the N-body case). 
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While the set iii) describes an infinite set of canonical relative variables with respect to the 
relativistic collective variables of the sets i) and ii), the sets i) and ii) describe a monopole 
field configuration, which depends only on 8 degrees of freedom like a scalar particle at 

rest [Pff, ~ 0] with mass M s = 



'(PJ) 2 — P| pa P^ but without the mass-shell condition 
M s pa const., corresponding to the decoupled collective variables of the field configuration. 
The conditions H(r, g) — K.(r,q) =0 select the class of field configurations, solutions of the 
Klein-Gordon equation, which are of the monopole type on the Wigner hyperplanes 



<pmon(r,a) = 2 J dq^F-(q)cu(q)P; - F(q)q- P+cos [q- (a - X*) - u(q)(r - XV) 
^ 2 X [PT ( dq^F-{q)u{q)cos \q- (a - X+) - oo{q){r - XV) 



7r m on(r,a) = -2 J dq^F-(q)u(q)PT-F(q)q-P^n^q-(a-X <j} )-Lu(q)(r-X T 4> ) 
pa -<lJp~l [ dqy/F-(q)u(q)sin \q ■ (a - X+) - u{q){r - XV) 



(3.69) 



If we add the gauge-fixings X^ pa to P<j, pa (this implies A(r) = and Hp = M^) and 
go to Dirac brackets, the rest-frame instant-form canonical variables of the Klein-Gordon 
field in the gauge r = T s are (in the following formulas it holds T s — X^ = — XV) 



a(T s ,q) = ^Fr(q)u;(q)P; + V ? li(T s ,q) 

e i[a{q)Xl+q-a}+i J dk f dk' K(T s ,k)g(k,k') A(k',q) 



iV = 5— + / dqV ? H(T s ,q), 
m J 



ir(T s ,a) 



<f>(T s ,a) = J dq^F-(q)u(q)P; + V € H(T 8 , q) 

' e i^(q)X^+q-a]+iJdkJdk'K(T s ,k)g(k,k')A(k',q) + 
_|_ e -i[u,(q)XT+q-a]-if dk J dk' K(T 3 ,k)g(k,k')A(k' ,q) 



= 2 J dq A 9 -(T S ; PT, H] cos \q-a + B,-(T S ; XJ, K 
= -i J dquj(q)^F-(q)uj(q)P; + V^Si(T 8 , q) 

' e i[oj(q)X^+q-a]+i J dk J dk' K(T S ,k)Q (k,k') A (£',?) _ 
_ e -i[^( (? )XJ+g-5 ! ]-i/c(fc/*'K(T s ,fc)g(fc,fc')A(fc',g)j = 

= -2 J dqu{q)A^T s - PT, H] sin [q- d + B^T 3 ; % K] ] , 
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A q iT s ;P;,B] = y/FT(q)u(q)P} + V?H(T s ,q) = y/I^j, 

Bf(T s ;X;,K] = J dkdk'K(T s ,k)g(k,k)A(k\fi+u(q)X; = 

= ^{T s ,q)-uo{q)T s . (3.70) 

The Hamiltonian Hp = M<p = P^ generates the following evolution in T s (= means 
evaluated on the equations of motion) 

° XI = {XL PI} = -1, =► XI = - T s , 



d 



dT f 



s 



° H(T„g) = 0, ^K(r si5 -) = 0, 



dT s v ' dT s 

J A 3 -(T S ; PI, H] = JrB^T.; XJ, K] = 0, 



<9T S yv " J <9T S 



d d 

>(T S , it) = - — — 0(T S , (?) = ir(T s , a), 



dT s Ty ' ' dX v 

d d 

?r(T s , <t) = - — — tt(T s , a) = -[A + m 2 ]0(T s , 5), 



d 2 d 2 

' +m 2 )0(T s ,a)=O. (3.71) 



dT? da 2 



- s 



Therefore, in the free case H(T s ,q), K(T s ,g*) are constants of the motion (complete 
integrability and Liouville theorem for the free Klein-Gordon field). Since the canon- 
ical variable P£ is the Hamiltonian for the evolution in T s = r, we need the inter- 
nal variable XJ = r + XT, (i.e. the internal time variable XTj to write Hamilton's 
equations -^rF = {F, P^} = —-§^r = — in the free case we have ^fr = — -^r 
on <f>(T s , a)[X^, P£, H, K] and tt(T s , <t)[XJ, P£, H, K], so that the evolution in the time 
XT = T s + XT, which takes place inside the Wigner hyperplane and which can be in- 
terpreted as an evolution in the internal time XJ, is equal and opposite to the evolution in 
the rest-frame time T s from a Wigner hyperplane to the next one in the free case. 

By adding the two second-class constraints X£ — T s = XT, ps 0, P£ — const. ~ 0, and by 
going to Dirac brackets, we get the rest-frame Hamilton- Jacobi formulation corresponding 
to the given constant value of the total ^energy: the field <f)(T s ,a), which is T s -independent 
depending only upon the internal time XT, becomes now even XT- independent. We find in 
this way a symplectic subspace (spanned by the canonical variables H, K) of each constant 
energy (P£ = const.) surface of the Klein-Gordon field. Each constant energy surface is 

not a symplectic manifold. However, it turns out to be the disjoint union (over XTj of the 

symplectic manifolds determined by XT = const., PI = const. 
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While the definitions of Subsection IIIC of the three external 4- and 3-center of mass is 

— * — # 

not changed, the three internal 3-centers of mass q^, R<f, and of Subsection HID have to be 
built by using the internal Poincare' generators (3.64) in Eqs. (3.16)-(3.19). In particular, 

the boosts are = X^P^ - M^X^ + Ks,/— - M^R^ with K r S lt> = -S™. Again, P^ pa 

implies ~ ~ Y$, the gauge fixing q^ pa implies pa 0, A(r) = 0, and the selection 
of the centroid (3.20) that coincides with the external Fokker-Pryce 4-center of inertia. 
The 4-momentum of the field configuration is peaked on this world-line while the canonical 
variables H (r, q), K (r, q) characterize the relative motions with respect to the monopole 
configuration of Eq.(3.69) describing the center of mass of the field configuration. 

As in the N-body case (ff + ^ q + ), the canonical 3-center of mass does not coincide with 
Xfj,, which in turn could be better defined as the center of phase of the field configuration. 
While the gauge fixing X^ pa 0, used to get Eqs. (3.70), implies pa Kg^ ^ 0, the gauge 

— * — * 

fixing pa implies X^ pa —Ks^/M^ (in fact it is Miller's definition of the 3-center of 
energy but only in terms of the spin part of the boost). 

Note that, like in Subsection HIE for the N-body case, there should exist a canonical 
transformation from the canonical basis = P£, X^, P^ pa 0, H(r, k), K(r, k), to a 

new basis gj, = yj (P^) 2 - P| pa (since {$/,, P£} = P<p/P^ pa is only weakly zero), 

q<j>, P<p ~ 0, H 9 (r, k), K 9 (r, k) containing relative variables with respect to the true center 
of mass of the field configuration 



Xtf, 


— * 

X<f> 


H(r, k) 




P<t>~0 


K(r, k) 



95 




H q (T,k)l 


Ml Pa M + 




K q (r,k)l 



(3.72) 



It does not seem easy, however, to characterize this final canonical basis. Besides the 
extension of the Gartenhaus-Schwartz methods from particles to fields (so that H q ~ H, 
K 9 ~ K), the real problem is finding the new internal time variable q^. 

In final canonical basis we would have still Eqs. (3.70) but with H g (r, k) and K ? (r, k) 
replacing H(r, k) and K(r, k) and with q^ replacing X^ as internal time. Since both the 

— * — * 

gauge fixings q^ ~ and X^ pa give A(r) = 0, both of them identify the same centroid 
(3.20) but lead to different internal times and relative variables connected by the canonical 
transformation (3.72). 

It is under investigation [52] the problem of characterizing the configurations of the real 
Klein-Gordon field in terms of dynamical body frames and canonical spin bases, to the 
effect of finding its orient at ion- shape variables. Ref. [26] contains also the treatment of the 
coupling of the real Klein-Gordon field to scalar particles and that of the charged (complex) 
Klein-Gordon field to electro-magnetic field. The collective and relative variables of the 
electro-magnetic field are now under investigation: they could be relevant for the problem of 
phases in optics and laser physics [53]. Finally, Ref. [27] contains the analysis of relativistic 
perfect fluids along these lines. 
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IV. THE MULTIPOLAR EXPANSION. 



In practice one is neither able to follow the motion of the particles of an open cluster in 
interaction with the environment nor to describe extended continuous bodies (for instance a 
satellite or a star), unless either the cluster or the body is approximated with a multi-polar 
expansion (often a pole-dipole approximation is enough). 

In this Section, after the treatment of the non-relativistic case in Subsection A, we show 
that the rest-frame instant form provides the natural framework for studying relativistic 
multipolar expansions of N free particles (Subsection B), of open clusters of particles (Sub- 
section C) and of the classical Klein-Gordon field (Subsection D) as a prototype of extended 
continuous systems (perfect fluids, electro-magnetic field,...). 



A. The non-relativistic case. 

The review paper of Ref. [29] contains the Newtonian multipolar expansions for a con- 
tinuum isentropic distribution of matter characterized by a mass density pit, a), a velocity 
field U r (t,<r), and a stress tensor <J rs (t, a), with p(t,a)U{t,a) the momentum density. In 
case the system is isolated, the only dynamical equations are the mass conservation and the 
continuity equations of motion 



dp(t,ff) dp(t,a)U r (t,a) 

dt da r 
dp{t,a)U r (t,a) _ d[pU r U s -a rs }{t,a) 

dt dV s 



0, (4.1) 



respectively. 

This description can be adapted to an isolated system of N particles in the following way. 
The mass density 

N 

p(t,a) = J2m t 5 3 (a-ff l (t)), (4-2) 
i=i 

satisfies 

= -f2mMt)-d^a-m)= f ^[pUl(t,al (4.3) 
i=i 

while the momentum density (this can be taken as the definitory equation for the velocity 
field) is 

N 

p(t,a)U r (t,a) = J2m l ff l (t)S 3 (ff ~m)- (4-4) 
i=i 

The associated constant of motion is the total mass m = YliLi- 

Introducing a function ((a, ffi) concentrated in the N points ffi, i=l,..,N, such that 
((a, ffi) = for a ^ ffi and ((ffi,ffj) = is a limit concept deriving from the char- 

acteristic function of a manifold), the velocity field associated to N particles becomes 



47 



N 

u( t ,a) = j2ma*,m)- (4.5) 
i=i 

The continuity equations of motion are replaced by 

8 f) N N 

-w, a)u r (t, a)} = -^- s Y. m rtm{tmt - m) + E = 

i=l i=l 

def d[pU r U s -a rs ](t,a) 

= fa • (46) 

For a system of free particles we have ff^t) =0 so that a rs (t, a) = 0. If there are inter- 
particle interactions, they will determine the effective stress tensor. 

Consider now an arbitrary point fj(t). The multipole moments of the mass density p, the 
momentum density pU and the stress-like density pU r U s , with respect to the point ff(t), are 
defined by setting (N > 0) 



m ri - r 



N 

N 

n = m[f)(t)] —m — ^mj, 

i=i 

p ri - r " r [fj(t)] = J d 3 a[a ri - 77 ri (t)]...[<7 r " - r) rn (t)]p(t,a)U r (t,a) = 

N 

= E^^MHO - v ri (t)]..M n (t) - v r -(t)], 
i=i 

N N 

n = P r [v(t)] = E m Mt) = E K i = K + ~ °> 

i=l i=l 

= E^w^w^rw -^(oi- (4-7) 

i=l 

The mass monopole is the conserved mass, while the momentum monopole is the total 
3-momentum, vanishing in the rest frame. 

If the mass dipole vanishes, the point fj(t) is the center of mass: 



N 

i=i 
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The time derivative of the mass dipole is 

= p r [ff(t)} - mfj r (t) =k t + - mfj r (t). (4.9) 

When ff(t) = q nr , from the vanishing of this time derivative we get the momentum-velocity 
relation for the center of mass 

p r [<lnr] = = ™?+ [~ in the rest frame}. (4-10) 
The mass quadrupole is 

N 

m rs [f](t)} =J2 m iV r i(t)v!(t)-rn V r {t)ri s {t) - (rf \t)m s [ff(t)} + V s (t)m r [rj(t)]) , (4.11) 
i=i 

so that the barycentric mass quadrupole and tensor of inertia are, respectively 

N 



m rs [q nr ] = ^2rriiriS(t)rit(t) -mq r nr q s nr , 
i=i 

I rS [Qnr] = 5 rS J2™ UU [qnr]-m rS [q nr } = 

u 

i=i 

1..JV-1 



a,b 

N-l 

m rs [q nr ] = 5 rs £ k <* Pa ■ Pb ~ I rS [Qnr] ■ (4. 12) 

a,b=l 

The antisymmetric part of the barycentric momentum dipole gives rise to the spin vector 
in the following way 

N N 

i=l i=l 
AT-1 

= ^Vl&r] = X ( 413 ) 

The multipolar expansions of the mass and momentum densities around the point 7/(t) 



are 



?1=0 

p(t, *)£T(f, a) = £ P w , ^ " V(t)). (4.14) 

n=0 
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Finally, we get the barycentric multipolar expansions as 



p(t,a) = m^(a-g nr )- V*[g^ 



da r da s 



n=3 



m ri ~~ rn [q nr \ d n 

n\ da r K..da r " 



5 3 (a-q nr ), 



p(t,a)U r (t,a) = K l5 3 (a-q nr )+[-e rsu S u +p^[q n 



_d_ 

do 7 ' 



5 3 {a -q nr ) + 



n....r n r r 

V { l>>r\ 



Qr, 



n=2 



n\ da r K..da r " 



S 3 (a-q n 



(4.15) 



B. The relativistic case on a Wigner hyper-plane. 



It is shown in Ref. [12] that, on a Wigner hyperplane with r = T s , the energy-momentum 
of N free scalar particles has the form 

T'"[ x P(T.)+4(u( Pa ))<T»] = e^u( Ps ))eUu(Ps))T AB (T s ,a) = 

N 



+ kr(T s ){u»(p s y r (u(Ps)) +u l/ (p s )e?(u(p s ))) + 



+ 



e»{u{Ps)KHp s )) 



N 



i=i 

N 

T Tr (T s , a) = J2 s3 (°-Vi(T s )K(T s ), 



i=i 

N 



T rs (T s ,a) = J2S 3 (°-UT S )) 



P» = tf = M sys u»(p s ) + €?(u(p a ))K r + « M S2/s ^(p s ), 

" I 



(4.16) 



We will now define the special relativistic Dixon multi-poles on the Wigner hyperplane 
with T s — t = for the N-body problem (see the bibliography of Ref. [11] for older attempts) . 
Note that, since we have not yet added the gauge fixings q + ~ 0, the centroid origin of the 3- 
coordinates has the form x^(r) = £s^ + ^(t) + J* r dri A r (rx) according to Eqs.(3.7) and (3.20). 
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Consider an arbitrary time-like world-line u> m (t) = z m (t, fftr)) = x^(t) -\-e%{u{p s )) r/ r (r) 



Xs" T '"(r) + equips)) fj r (r) [fj r (r) = ^ r (r) + J q t c?ti A r (ri)] and evaluate the Dixon multi-poles 
[28] on the Wigner hyper-planes in the natural gauge with respect to the given world-line. 
A generic point will be parametrized by 



z»(r,a) = x^T) + ^(u( Ps ))a r = 

= w^{t) + e?{u{p s ))[a r - V r (r)] =V(r) + 8z»(t, a), (4.17) 

so that 5z^(t, a)u^(p s ) = 0. 

While for fj{r) = [rj(r) = dr\ \ r {j\)\ we get the multi-poles relative to the centroid 
rr^(r), for r/(r) = we get those relative to the centroid x^ + ^(r). In the gauge R + ~ q + ~ 

— * — # 

y + pa 0, where A(r) = 0, it follows that fj{r) = r/(r) = identifies the barycentric multi-poles 
with respect to the centroid xF + ^ M (t), that now carries the internal 3-center of mass. 

Lorentz covariant Dixon's multi-poles and their Wigner covariant counterparts on the 
Wigner hyper-planes are then defined as 



6^(^(p s ))...6^(t i ( Ps ))6^( W (p s ))6^( W (p s ))^-- AB (7;, ? 7) 



d 3 a5z^(T s ,a)...5z^(T s ,a)T^[xf^(T s ) + e ^u{p s ))a u } = 
= e» A (u(p s )y B (u( Ps )) J d 3 a5z^(T s ,a)....5z^(T s ,a)T AB (T s ,a) = 

= e£(«(p.))-#(«(p.)) 

r N l 



i=l 



+ 6?(u(p s )y s (u( Ps )) 



N 



N 



n...r n AB 



(T S ,V) = J d 3 a[a r ^ - V ^(T s )]...[a r - - V r -]T AB (T s ,a) = 
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N 

= W EfoPCT.) - V n (T.)]...[v?(T.) ~ V rn (T s )]y/m> + %(T.) + 



i=i 

AT 



i=i V" l i+'W«J 



A? 



+ (« + W) " ^m)]...k r "(T s ) - 77- (T S )]<(T S ), 



1=1 



u,M ti?~^(T s ,ff) = 0. (4.18) 

Related multi-poles are p£-^(T s ,rf) t^-^ u (T s ,rf)u u (p s ) 
^(u(p s ))..<(«(p a ))^( U (p a ))g?"- r » AT (7;,J7). They satisfy u^{p s )p^-^{T s ,ff) = 
and for n = they imply p^(T s , if) = e^(u(p s ))q^ T (T s ) = P£ pa p^. 

The inverse formulas, giving the multipolar expansion, are 



T^K(T S ) + 6zP(T s , a)} = T^[xf^(T s ) + e?(u(p s )) a r \ = 

= e\Hp s )y B (u(p s ))T AB (T Sj a) = 

oo ri...r„AB/rp an 



n=0 

Note however that, as pointed out by Dixon [28], the distributional equation (4.19) is 
valid only if analytic test functions are used defined on the support of the energy-momentum 
tensor. 

The quantities q^- rnTT (T s ,ff), q^- rnrT (T s ,ff) = q^~ rnTr {T s ,ff), q^- rnUV {T s ,ff) are the 
mass density, momentum density and stress tensor multi-poles with respect to the world-line 
w^{T s ) (barycentric for ff = fj — 0). 



1. Monopoles 



The monopoles, corresponding to n = 0, have the following expression (they are rf- 
independent) 



N 



q AB (T s ,ff) = 5*5* M + 5 A 5 B —[ 



+ (s A S B + 5 A 5 B 



mf + k; 



)«: 



u 

+' 
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N 



q T T T (T s ,ff) 

q r T T (T s ,v) 
Qta( T s,v) 



J2 m iC 2 + H re i + 0(l/c), 



i=i 



= k T j^ ~ 0, rest — frame condition (also at the non — relativistic level) , 



1..N-1 



AT 9 

mf 



N 

■ c _oo ^rriiC 2 - H rel + 0(1/ c). 



(4.20) 



i=i 



Therefore, independently of the choice of the world-line w^(t), in the rest-frame instant 
form the mass monopole q^ T is the invariant mass M sys = J2iLi \/ m i + while the mo- 
menta monopole q r ^ vanishes and g™ is the stress tensor monopole. 



2. Dipoles 



The mass, momentum and stress tensor dipoles, corresponding to n — 1, are 



+ 



q r T AB (T s ,v) = S^M[Rl(T s )-n r (T s )]+5X\j2 ^M T *) ~ V r (T 8 )q^(T a ,ff) 

N 



+ 



i=l 



The vanishing of the mass dipole qJ r rT implies f/(r) = n(r) — J* dri A(ti) = R + and 
identifies the world-line #(r) = xF + ^(r) + e(f(u(p s )) + f r dn A r (ri) . In the gauge 

i? + m q + y + m 0, where A(r) = 0, this is the world-line w^(t) = x^ + ^(t) of the 
centroid associated with the rest-frame internal 3-center of mass q + . We have, therefore, 

— * 

the implications following from the vanishing of the barycentric (i.e. A(r) = 0) mass dipole 



q r T ^(T s , ff) = e£ (u(p s ))t£ (T s , if) — M [R r + (T S ) - n r (T s )] = 0, and A(r) = 0, 



ff(T s ) = jj(T 9 ) =R + ^ q + » y + . 



(4.22) 
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In the gauge R + pa q + pa y + pa 0, Eq.(4.22) with fj — fj — implies the vanishing of the 
time derivative of the barycentric mass dipole: this identifies the center- of -mass momentum- 
velocity relation (or constitutive equation) for the system 



(11 e 



(4.23) 



The expression of the barycentric dipoles in terms of the internal relative variables, when 
fj = fj = R + pa q + pa and k + pa 0, is obtained by using the Gartenhaus-Schwartz transfor- 
mation. 



q r T TT (T s ,R + ) = 0, 



N 



N-1 



q r T UT (T s ,R + ) 



E « - r>\ = E p>« + « - 

i=l a=l 
N-l 1..N-1 

E«= E 



q r T uv (T s ,R, 



a=l 

N 

E^ 



«6 



K u v 



JV 



i=l 



ft 



-*;E 



V JV-1 



1=1 



ft 



K U K V 



=1 a=l * i=l " l 

.JV-1 JV v^AT 

Er AT lailbilci 2^j=i m j1ajl 
rE— — jwv+o(i/ C ). 



abc i=l 



(4.24) 



The antisymmetric part of the related dipole pf^ lfl {T s , fj) identifies the spin tensor. Indeed, 
the spin dipole is 



S^(T s )[fj\ = 2p [ r\T s ,ff) = 2^(u(p s ))ei(u(p s ))q^(T s ,fj) = 

= M sys [R r + (T S ) - 7f (T s )] \e»{u(p s ))u v {p s ) - <(u(p 8 )K(p a ) 



+ 



N 



i=l 



+ E^w^kto c{f(«(p.))<(«(p.)) - <(«(p.))c?(«(p.)) 



K(p s )( T s,V) = u,(p s )S^(T s )[fj\ = -e:(u( Ps ))[S: r -M sysV r (T s )] = 

= -e»{u{Ps))M sys [R r + (T S ) - V r (T s )} = -e»(u(p s ))q r T ^(T s ,fj), 



u,(p s )S^(T s )[fj\ = 0, ^v = R + , 
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Jj bary centric spin for k + k 0, ff—fj — 0, see Eq(4.17), 
SrCW=0] = S^ = e rsu S:e?(u(p s )K(u(p s )). (4.25) 



This explains why m^ p ^(T s ,ff) is also called the mass dipole moment. 

We find, therefore, that in the gauge R + w <f + p» y + 0, with ss M sys u fM (p s ) = 
M sys x^s +y>ll (T s ), the M0ller and barycentric centroid xi' + ^(T a ) is simultaneously the Tul- 
czyjew centroid [48] (defined by S**" P u = 0) and also the Pirani centroid [47] (defined by 
xit + ^ = 0). In general, lacking a relation between J^-momentum and ^-velocity, they are 
different centroids. 



3. Quadrupoles and the barycentric tensor of inertia 



The quadrupoles, corresponding to n — 2, are 



N i 

q r T ^ AB (T s , ff) = 5 A 5 B 5>p(T.) - V ri (T s )M 2 (T s ) - ^{T s )\ yjm* + £f(T s ) + 

i=i 

JV 

+ « Ek n m) - ^ ri m)]k-(T s ) - ^{t s )\ ^pL= m) + 

JV 

+ (« + W) E^m) - ?f H^MTO - ?f 2 (T s )K(T,), 
i=i 

(4.26) 

When the mass dipole vanishes, i.e. ff = R + = ^ ffi \[m\ + K%/M sys , we get 



JV 



i=l 
N 



i=i 



JV 



q r ™ v (T s ,R + ) = -WW- Z?) 



i=i 



V™1 + Ws) 



1..JV 1..N-1 



— E ^ a ^' 



(4.27) 
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Following the non-relativistic pattern, Dixon starts from the mass quadrupole 

N l 

q^ TT (T s , R + ) = Et« 2 V m * + %](T S ) - M sys i^, (4.28) 
i=\ 

and defines the following barycentric tensor of inertia 

QZniTs) = ^ E Qt UTT (T s , R + ) - q r T ^ TT (T s , R + ) = 



u 

N 



1..N-1 

-c^oo E k ^ ■ P^ riT2 ~ P>2] + (VC) = r ir2 [qnr] + 0(1/C). 
ab 

(4.29) 

Note that in the non-relativistic limit we recover the tensor of inertia of Eqs.(4.11). 
On the other hand, Thome's definition of barycentric tensor of inertia [30] is 

QlUTs) = *™ J2 Vt A a(T s , R + ) - q^ A A (T s , R + ) = 

u 

j2 m K^ r2 (v t - R + f - (y? - R7M 2 - R7)) {Tg) 



E k *\P«> ■ P^ rir2 - P> T $\ + °(VC) = I Tir2 [qnr] + 0(1/C). 



i=l 
1..JV-1 



c^oo 

(4.30) 

In this case too we recover the tensor of inertia of Eq.(4.11). Note that the Dixon and 
Thorne barycentric tensors of inertia differ at the post-Newtonian level 



^ 1..JV-1 l.JV m m 

^dixoni^s) ~ hhornei^s) = ~ (7ai — 7aj ) (7f« ~ 7f>fc) 

ab ijk 



\p qa ■ P qb S r ^ - P r qaP r q l\ lcadl7lqc ' ^ + 0(l/c 2 ). 

(4.31) 
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4- The multipolar expansion 



It can be shown that the multipolar expansion can be rearranged in the form 



TV[ x if + )0(T.) + e?(u(p.)K] = T^[w%T s ) + ^ r {u{Ps)W ~ V r (T s ))} = 
= u^(p s y2(u(p s ))[5?M sys + tf<tf S (* " V(T S )) + 
+ \s^{TMu v \Ps)e^u{p s ))^- r 5\a - ff(T s )) + 

00 ( — ~\\ n P) n 

+ E 4r^'^ r " ^ («(P.))-^("(P.)) ^ ^ *V " ^)), (4.32) 

n=2 

where for n > 2 and 77 = 0, we have I^-^ U (T S ) = £M j^«-il^->(T s ), with 
j^i.-nnnvpa being the generalized Dixon 2 2+n -pole inertial moment tensors given in 
Ref.[ll]. 

Note that, for an isolated system described by the multi-poles appearing in Eq.(4.32) 
[this is not true for those in Eq.(4.19)] the equations d^T^ — O imply no more than 
the following Papapetrou- Dixon- Souriau equations of motion [29] for the total momentum 
P%(T S ) = ^ A {u(p s ))q^ T (T s ) = p» and the spin tensor S^(T s )[ff = 0] 



dP£(T s ) 
dT x 



= 0, 



dSP(T )[ff-0] o 2p \P {TsK]{ps) = 2 ^{u{p s ))u»\p s ) » 0, 



dT s 



^ ^=0, ^±=0, ^0. (4.33) 
dl s dl s dl s 

5. Cartesian Tensors 

In the applications to gravitational radiation, irreducible symmetric trace-free Cartesian 
tensors (STF tensors) [30, 31] are needed instead of Cartesian tensors. While a Cartesian 
multi-pole tensor of rank I (like the rest-frame Dixon multi-poles) on R 3 has 3' components, 
|(/ + 1)(/ + 2) of which are in general independent, a spherical multi-pole moment of order I 
has only 2/ + 1 independent components. Even if spherical multi-pole moments are preferred 
in calculations of molecular interactions, spherical harmonics have various disadvantages in 
numerical calculations: for analytical and numerical calculations Cartesian moments are 
often more convenient (see for instance Ref. [54] for the case of the electrostatic potential). 
It is therefore preferable using the irreducible Cartesian STF tensors [55] (having 2/ + 1 
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independent components if of rank I), which are obtained by using Cartesian spherical (or 
solid) harmonic tensors in place of spherical harmonics. 

Given an Euclidean tensor A kl ... kl on R 3 , one defines the completely symmetrized tensor 
S kl __ kl = A( fcl .. fej ) = j\'52 ir A k , 1) _ k , I y Then, the associated STF tensor is obtained by 
removing all traces ([1/2] = largest integer < 1/2) 



[J/2] 



+ l...k I )irh..-jnjn1 



n=0 



-l) n 



l\(2l-2n-l)\\ 



For instance (T abc ) STF = T {abc) - \ 



(l-2n)\(2l-l)\\(2n)\\ 

0~abT(iic) + 0~acT(ibi) + h C T( a ii) 



(4.34) 



C. Open N-body systems. 

Consider now an open sub-system of the isolated system of N charged positive-energy 
particles plus the electro-magnetic field in the radiation gauge (see the second paper of 
Ref. [6]). The energy-momentum tensor and the Hamilton equations on the Wigner hyper- 
plane are, respectively, [to avoid degenerations we assume that all the masses are different; 

7?_l = EA 



N I 

T^(r,a) = ^\^-Ur))\lrnj+[K l (T)-Q l A ± (r,m(rW + 



N d 1 

+ ^Qi s ) x — 5 3 (a - f)i(T)) + - [ttI + B 2 ](t, 3) + 

i=i 



i,k,i=ik 



N 



r-(r,a) = Y,^^-V l (r)M(r)-Q t A r L (r,ff t (r))] + 



i=l 



+ [{^± + J2Q^ S3 ^ -Vi(r))) x B](r,a), 
i=i 



m?+ [«i(r) -QiA ± (r,77i(r))] 2 



?=i 
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N f) N f) 

- [(*!_ + £ ft A 53(C? _ ( ?± + S ft " * (t)) )' 

i=l i=l 

+ S r S a ]](r,5). (4.35) 



Kj( T ) ~ QiA±(T,f)j(T)) 



, . ^ Q^(r) - y r)) + g? a ^ 

^ 4vr | 77i(r) - 77*. (t) | d c% 



nl(T,a) = AA r ± (r,a)-^QiP r ±^mr)S 3 ^-Vi(r)l 

i 

k + (t) + J d 3 a[rr± x B](t, a) ~ (rest — frame condition). (4.36) 

Let us note that in this reduced phase space there are only either particle- 
field interactions or action-at-a-distance 2-body interactions. The particle world-lines 
are xf(r) = + u^(p s )r + e%(u(p s )) ?7[(t), while their 4-momenta are pf(r) = 

^/ m 2 + ^._g.^ ±(r ^.)]2^( ps ) + <£( u (p a )) [«r - A r ± (r, ff t )}. 
The generators of the internal Poincare group are 

N I 

V\ mt) = M = J2\Jrf + (Ki(r) - QiA ± {T,Tk(T)))* + 

i=i 

P(mt) = «+(r) + y dV[7T_L x J3] (r, a) « 0, 

^ r 

^(int) = J2^i( T ) x Ki(r)) r + / d 3 a(ax [7r ± xfi] r (r, a), 
i=i ^ 

* I 

fc(int) = - ^(r)^m^+[K t (r)-Q t A ± (r, ^r))f+ 



i=i 
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N 1..N 



+ \ E E °-> I d 3 aa r c(a-ff l (T))-c(a-ff J (r)) + 

i=l i+i J 

+ Q t J d 3 a7T r ± (r,a)c(a- ^(r))] - \ J d 3 aa r + B 2 )(r, a), (4.37) 

with ciffi-ffj) = 1/(4^1^-^1) [Ac(<t) = 5 3 (a), A = -d 2 , c(a) = dc(a) = a/(4vr |a| 3 )]. 
^(int) = 1 TT an d T^lint) = Q TT are the mass and momentum monopoles, respectively. 

For the sake of simplicity, consider the sub-system formed by the two particles of mass m x 
and rri2- Our considerations may be extended to any cluster of particles. This sub-system is 
open: besides their mutual interaction the two particles have Coulomb interaction with the 
other N — 2 particles and are affected by the transverse electric and magnetic fields. 

Exploiting the multi-poles we will select a set of effective parameters (mass, 3-center of 
motion, 3-momentum, spin) describing the two-particle cluster as a global entity subject to 
external forces in the global rest-frame instant form. This was indeed the original motivation 
of the multipolar expansion in general relativity: replacing an extended object (an open 
system due to the presence of the gravitational field) by a set of multi-poles concentrated on 
a center of motion. Now, in the rest-frame instant form it is possible to show that there is 
no preferred centroid for open system so that, unlike the case of isolated systems where, in 
the rest frame k+ 0, all possible conventions identify the same centroid, different centers 
of motion can be selected according to different conventions. We will see, however, that one 
specific choice exists showing preferable properties. 

Given the energy- momentum tensor T (r, a) (4.35) of the isolated system, it would seem 
natural to define the energy-momentum tensor T^(r, a) of an open sub-system composed 
by a cluster of n < N particles as the sum of all the terms in Eq.(4.35) containing a 
dependence on the variables fji, K iy of the particles of the cluster. Besides kinetic terms, 
this tensor would contain internal mutual interactions as well as external interactions of 
the cluster particles with the environment composed by the other N — n particles and by 
the transverse electro-magnetic field. There is an ambiguity, however. While there is no 
problem in attributing to the cluster the whole interaction with the electro-magnetic field, 
why should we attribute to it just all the external interactions with the other N — n particles! 
Since we have 2-body interactions, it seems more reasonable to attribute only half of these 
external interactions to the cluster and consider the other half as a property of the remaining 
N — n particles. Let us remark that considering, e.g., two clusters composed by two non- 
overlapping sets of ni and n2 particles, respectively, since the mutual Coulomb interactions 
between the clusters are present in both and T^?), according to the first choice we 

would get T c ^ ? f i+n2 - ) 7^ + T^?y On the other hand, according to the second choice 

we get T^ i+n ^ = + T^ 2 y Since this property is important for studying the mutual 

relative motion of two clusters in actual cases, we will adopt the convention that the energy- 
momentum tensor of a n particle cluster contains only half of the external interaction with 
the other N — n particles. 

Let us remark that, in the case of /c-body forces, this convention should be replaced 
by the following rule: i) for each particle rrii of the cluster and each /c-body term in the 
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energy-momentum tensor involving this particle, k = hi + (k — hi), where hi is the number 
of particles of the cluster participating to this particular /c-body interaction; ii) only the 
fraction hi/k of this particular /c-body interaction term containing mi must be attributed to 
the cluster. 

Let us consider the cluster composed by the two particles with mass m\ and 777.2- 
The knowledge of T^ B d = T^S on the Wigner hyper-plane of the global rest-frame instant 
form allows us to find the following 10 non conserved charges [due to Ql — we have 

Jmj + [Ki - Qi A ± (r, th)] 2 = + % ~ Qi ^=#1 



2 

f d 3 aTr(r,a) = £ \J mf + [^(r) - Q t A±(r, f]i(r))] 2 + 

"* i=i 

2 

Q1Q2 1 <5i Qfc 

47r|r7 1 (r)-r7 2 (r)| 2 2 ^ ^ 4tt ^(r) - r/ fc (r)| 2 " 

= M^t) + M c ( ext ), 

II ,f / 2 1 -2 Q1Q2 

^ V 471(77! (r) -?72(t)| 2 

Vc = {J d 3 aT r c T {r, a)} = k ± (t) + k 2 (t), 

Jc = {e ruv J d 3 a[a u T^ -a v T: T ](r,a)} = ff l (r)xK 1 (r)+ff2(r)y<Mr), 

K, c = - f d 3 aaTr(r,a) = - ]T r/,(r) \J m 2 + [^(r) - l ± (r, ^(r))] 2 - 
^ i=i 

~ y^.Qj I d 3 an ± (T,a)c(a -ffi(T))- 
i=i ^ 

- Q1Q2 J d 3 a a c(a - ff^r)) ■ c(a - fj 2 (r)) - 

1 2 /" 

- - Q fc / d 3 crac(a -ffi(T)) ■ c\a -ff k (T)) = 

i=l fe^l,2 ^ 

— * — * 

fc c (i n t) = Vi( r ) \J m ! + - Qi Q 2 y rf 3 crac(a-r7i(r)) •c(a-r7 2 (r)). 

(4.38) 



i=l 



Such charges do not satisfy the algebra of an internal Poincare' group just because of the 
openness of the system. Working in an instant form of dynamics, only the cluster internal 
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energy and boosts depend on the (internal and external) interactions. Again, M c = (fj and 
V r c = q r J are the mass and momentum monopoles of the cluster. 

Another quantity to be considered is the momentum dipole 



P? 



2 2 . 

i£(t)k?(t) - E Qi / dVc(*-#(T))[FAi + 0Mi](T,*), 
1=1 1=1 J 

2 

i=i 

2^Q, / ^c(a-^(r))[^Al + 9 s Al](r,a), 
i=i ^ 

ru Mr rati rrv 

Pc Pc ~ t <J C - 



(4.39) 



The time variation of the 10 charges (4.38) can be evaluated by using the equations of 
motion (4.36) 



dM c 



= 5>( 

i=i 



Ki(r) • n±(T,ffi(T)) 



+ 



fc^l,2 



V 7 ^- + k- v 7 "^ + ^ 47r|^(r) -?/ fc (r)| 3 7' 



= E^ 



¥1,2 



47r|7/i(r) -7ffc(i 



E^ 



x A ± (r,r/i(r)) +rfi(r) x 



9 1 



A ± (r,77i(r)) - 



47r |r7i(r) -rf fc (r)| 3 /' 



^ = Ql Q 2 / rfW 



+ c(<t-?7i(t)) 



«i(r) 



A/mf + k2(t) 
«2(r) 



9) c(a-rfi(r)) -c((T-772(t)) + 



v/m| + k|(t) 



•9 c(a-7f 2 (r)) 
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fc#l,2 7 



i=l 

+ c(CT-r/i(r)) 



v 7 ™! + 4( r ) 



dj c(a-f) k (T)) ) 



9) c{a-ffi(r)) 



c(a-fj k (T)) + 

(4.40) 



Note that, if we have two clusters of n x and n 2 particles respectively, our definition of 
cluster energy-momentum tensor implies 



M c ( ni+n2 ) = M c ( ni ) + M c („ 2 ), 
^ > c(m+n2) ^c(m) ~l~ ^Pc(n2)i 

fcc( ni +n 2 ) = /C c ( ni ) + /C c („ 2 ). (4-41) 



The main problem is now the determination of an effective center of motion Q( T ) 
with world-line u>£(t) = a;^ + u^(p s )r + e(f(«(p s )) Cc ( r ) i n the gauge T s = r, q + = 

— * 

R + = y + = of the isolated system. The unit 4-velocity of this center of motion is 

<(r) = <(r)/\/ 1 - C(r) with <(r) = u»{p s ) + e?(u(p s )) £(r). By using 8 z"{t,3) = 
e r( u (Ps)) (°" r_ C r ( r )) we can define the multipoles of the cluster with respect to the world-line 

<(r) ' 

q7- rnAB (r) = J d'a [*" - C(t)]..[o*" - C(r)] T c AB (r, a). (4.42) 

The mass and momentum monopoles, and the mass, momentum and spin dipoles are, 
respectively 



q? = M c , C = V r c , 

qrjr = _ K _ ^ C(r) = ^ {K{t) _ C(r))) g rur = ^ _ ^ ^ 



Sr = K(u( Ps )) u»{p s ) - e»(u(p s )) u»(p s )} qT + e?(«(p.)) <(«(p.)) (?i 
= [e?(u(p.)) u"{p s ) - e»(u(Ps)) «"(p.)] M c (R r c - C) + 
+ e?(u(p s )) e v u {uip s )) [e™ Jl - (C V u c - Q K)] , 



rur 
c 



urr 
< dc 



) = 



m 



C(p s ) 



-S^u u {p s ) = -e»{u{p s ))q: 



(4.43) 
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Then, consider the following possible definitions of effective centers of motion (clearly, 
many other possibilities exist) 

1) Center of energy as center of motion, ( c (e)(t) = R c (t), where R c {t) is a 3-center of 
energy for the cluster, built by means of the standard definition 



3 _ £c 

J^c — TT- 

M r 



(4.44) 



It is determined by the requirement that either the mass dipole vanishes, q r c TT = 0, or the 
mass dipole moment with respect to u^(p s ) vanishes, ^c( Ps ) = 0- 

The center of energy seems to be the only center of motion enjoying the simple compo- 
sition rule 



R 



c(m+n2) 



^c(ni) Rc(ni) + M c ( m ) R c ( n2 ) 



M, 



c(ni+n 2 ) 



(4.45) 



The constitutive relation between V c and R c {t), see Eq.(4.23), is 



dq r c 



dr 



-K, r c -M c R r c -M c R r c , 



4 



J VLV A/mf + k((t) > 

Mr) 



+ C\(7-?7i(T 



))■[( 



d) c(a-ff 2 (r)) 



Eft E ft -, ==n-9)c(a 



fji{r)) ■ c(a - if k (T)) + 



+ c{a - Ur)) ■ \( . ff! 2 , . • d) c(a - fj k (r))} ) 
LV A/mt + kZ(t) ' I' 



(4.46) 



/,From Eq.(4.25), it follows that the associated cluster spin tensor is 



= e?(u(p s )) <(u(p.)) e r ™ \r c - (R c x V c y 



(4.47) 



2) Pirani centroid Cc(p)(t) as center of motion. It is determined by the requirement 
that the mass dipole moment with respect to 4- velocity wjf(r) vanishes (it involves the 
anti-symmetric part of p^ r ) 



64 



m 



C(w c ) 



Sc" ™cu = 0, =>- Cc(P) ■ Cc(P) = Cc(P) ' -^c, 



4 



Cc(P)(r) 



M c -V c -( c(p) (t 



M c R c -R c -( c{p) (t)V c -( c(p) {t)xJ c . (4.48) 



(P)> 



Therefore this centroid is implicitly defined as the solution of these three coupled first 
order ordinary differential equations. 

— * 

3) Tulczyjew centroid Cc(t) (t) as center of motion. If we define the cluster 4- momentum 

P£ = M c u^(p s ) + V S C e£(u(p a )) [P 2 = M 2 c -V 2 c d = M% its definition is the requirement that 
the mass dipole moment with respect to P 6 M vanishes (it involves the anti-symmetric part of 

p7) 



m 



c(Pc) 



-S^P a/ = 0, ^V c -( c(T) =V c -R c , 



4 



Cc(T)(t) 



Ml - VI 



M 2 R C -V C -R C V C -V C x J c 



(4.49) 



Let us show that this centroid satisfies the free particle relation as constitutive relation 



V c = M c ( 



c(T)» 



i * / 

Qc(T) 



M r 



Ml - VI 



V Z C R C + V C -R C V C + V C x J c 

Sr = WHPs)) u»{p s ) - equips)) u»{p s )\ £fc + 
+e?(u( Ps )) <(u(p.)) \T c - (Cc(T) x P> 



(4.50) 



If we use Eq.(4.48) to find a Pirani centroid such that ( c = V C /M C) it turns out that the 
condition (4.48) becomes Eq.(4.49) and this implies Eq.(4.50). 

The equations of motion 



M c (r) U(r) = V c (t) - M c (r) ( c(T) (r), 



(4.51) 



contain both internal and external forces. In spite of the nice properties (4.50) and (4.51) of 
the Tulczyjew centroid, this effective center of motion fails to satisfy a simple composition 
property. The relation among the Tulczyjew centroids of clusters with ni, n<i and n\ + n<i 
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particles respectively is much more complicated of the composition (4.45) of the centers of 
energy. 

All the previous centroids coincide for an isolated system in the rest-frame instant form 
with V c = k + ~ in the gauge q + ~ R + ~ y + ~ 0. 

4) The Corinaldesi-Papapetrou centroid with respect to a time-like observer with ^-velocity 
V ^( T )> Ca(cp)( T ) as center of motion. 



Clearly these centroids are unrelated to the previous ones being dependent on the choice of 
an arbitrary observer. 

5) The Pryce center of spin or classical canonical Newton-Wigner centroid ( c {nw)- 
It defined as the solution of the differential equations implied by the requirement 
{C(w)>C(w)} = °> {C( N wy V c} = & rs - Let us remark that, being in an instant form 
of dynamics, we have {V^,V^} = also for an open system. 

The two effective centers of motion which appear to be more useful for applications are 
the center of energy ( c (e)( t ) and Tulczyjew's centroid Cc(t)(t), with Cc(B)( r ) preferred for 
the study of the mutual motion of clusters due to Eq.(4.45). 

Therefore, in the spirit of the multipolar expansion, our two-body cluster may be de- 
scribed by an effective non-conserved internal energy (or mass) M c (r), by the world-line 
w c = x o + ufl {Ps) r + e r( u (Ps)) (c(EorT)( T ) associated with the effective center of motion 

Cc(EorT)ij) and by the effective 3-momentum P c (t), with Cc(EorT)ij) and V c {r) forming a 
non-canonical basis for the collective variables of the cluster. A non-canonical effective spin 
for the cluster in the 1) and 3) cases is defined by 

a) case of the center of energy: 



m 




-S^v v = 0. 



(4.52) 



Sc(e)(t) 



J c (t)-R c (t)xV c (t), 



d(c(E){r) 
dr 



dj c (r) dR c (r) 



dr dr 



x V c {t) - R c (t) x 



dV c (r) 



dr 



b) case of the Tulczyjew centroid: 



£ (T )(t) = J c (r) - Ce(T)(r) x V c (t) = 

M c 2 (r) S c(e) (t) - V c {r) ■ J c {r) V c {r) 



<Kc(t)(t) dj c {r) r dV c {r) 



(4.53) 
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Since our cluster contains only two particles, this pole-dipole description concentrated 
on the world-line w£(t) is equivalent to the original description in terms of the canonical 
variables ffi(r), Ki(t) (all the higher multipoles are not independent quantities in this case). 

Finally, in Ref. [12] there is an attempt to replace the description of the two body system as 
an effective pole-dipole system with a description as an effective extended two-body system by 
introducing two non-canonical relative variables p c {EorT)ij), ^c{EorT)ij) with the following 
definitions 



def p 1 / 1 . x 

Vl = ic(EorT) + ^ Pc(EorT), Qc(EorT) = ^ KVl + V2), 

-> def 1 -, -> -> -> 

V2 = U(EorT) ~ -Pc(EorT), Pc(EorT) — T)i — 7) 2 , 

—* def 1 ^1 —> ^1 — » — * 

«1 = ^ lc + ^c(EorT), P C = K 1 + K 2 , 

_ de/ 1 ^ _ 1 _ v 

Jc = rji X Ki + rj 2 X K 2 = (c(EorT) X V c + Pc(EorT) X ^c(EorT), 



Sc(EorT) — Pc(EorT) X ^c(EorT)- 



(4.54) 



Even if suggested by a canonical transformation, it is noi a canonical transformation and it 
only exists because we are working in an instant form of dynamics in which both V c and 
J c do not depend on the interactions. Note that we know everything about this new basis 
except for the unit vector p c (EorT)/\Pc(EorT)\ and the momentum Ti c {EorT)- The relevant 
lacking information can be extracted from the symmetrized momentum dipole p r c u + p" r , 
which is a known effective quantity due to Eq.(4.40). However, strictly speaking, this type 
of attempt fails, because p r c u + p c Lr does not depend only on the cluster properties but also on 
the external electro-magnetic transverse vector potential at the particle positions, as shown 
by Eq.(4.39). Consequently, the spin frame, or equivalently the 3 Euler angles associated 
with the internal spin, depend upon the external fields. 



D. The Multipoles of the Real Klein-Gordon Field. 

In the rest-frame instant form we have the following expression [26] for the energy- 
momentum of the real Klein-Gordon field 



T^[x»(T s ) + eZ(u(Ps))o- u M = T»V s (T s )+eZ(u(p s ))a u }[X£,P^H : K} = 

= l -u»{p s )u»{p s )[n 2 + 0<P) 2 + m 2 2 ](T s ,a) + 

+ ^Mp s )yMPs))[-\Srs[* 2 ~ (<90) 2 - m 2 2 ] + 
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+ <9 r 0<9 s 0](T s ,a) - 

[u^{p s )e v MPs))+u v {p s )e^Ps))]WWs^) = 
= p[0, n}u»( Ps X(p s ) + P[0, tt] [rT ~ u»{p s )u»{p s )\ + 

p[M = i[7r 2 + (90) 2 + m 2 2 ], 

V[M = ^[7T 2 -^(90) 2 -m 2 2 ], 
<H<M = -ird r <f>e?(u(p s )), 

^ ess m,,T)[0] = e;(^))et(«(p s ))T^[<(T s ) + eC( W (p s )) ( 7l[0] = 
= W<pd s <P](T s , a) - - (90) 2 - m 2 2 ](T s , a), 

p$[</>] = y , d 8 (7r^[^(r.) + cS(«( P .))(7«][^K( P- ) = 
= p;/(p s ) + py r {u{p s )) « p».) ~ rf, 

M = = P;. (4.55) 

The stress tensor T^ ress (T s , a)[(p} of the Klein-Gordon field on the Wigner hyper-planes 
acquires a form reminiscent of the energy-momentum tensor of an ideal relativistic fluid as 
seen from a local observer at rest (Eckart decomposition; see Ref. [27] ): i) the constant 
normal u^{p s ) to the Wigner hyper-planes replaces the hydrodynamic velocity field of the 
fluid; ii) p[0, n](T s , a] is the energy density; iii) V[(p, n](T s , a) is the analogue of the pressure 
(sum of the thermodynamical pressure and of the non-equilibrium bulk stress or viscous 
pressure); iv) g M [0, tt](T s , <?) is the analogue of the heat flow; v) T™ stress [(f), 7r](T s , a) is the 
shear (or anisotropic) stress tensor. 

Given the Hamiltonian version of the energy- momentum tensor, Eqs.(4.18) allow us to 
find all Dixon multi-poles of the Klein-Gordon field [26] , with respect to a natural center of 
motion identified by the collective variable defined in Subsection IIIH. Finally, Eqs.(3.70), 
conjoined with the assumption that we can interchange the sums with the integrals, allow 
us to define another class of multi-poles [26] with respect to the centroid Xg(r), origin of 
the 3-coordinates a, in the gauge ~ 0, which, when the fields have a compact support 
in momentum space, form a closed algebra with a generalized Kronecker symbol, that could 
be quantized instead of the Fourier coefficients. 
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V. FINAL COMMENTS AND OPEN PROBLEMS. 



In this article we have shown how the traditional Jacobi technique based on the clustering 
of centers of mass for an N-body system, can be profitably replaced by a technique based 
on the clustering of spins in order to develop a theory of the relativistic, rotational, and 
multi-pole, kinematics for deformable systems. More generally, the relativistic extension is 
also made possible by a systematic use of the so-called rest-frame instant form of relativistic 
dynamics on Wigner hyper-planes, orthogonal to the conserved 4-momentum of the isolated 
system. This is a form of dynamics characterized by a typical doubling of the Poincare 
canonical realizations into so-called external and internal realizations. In fact, this framework 
appears to be the most natural theoretical background for the description of isolated systems 
(particles, strings, fields, fluids) in special relativity. 

The rest-frame instant form of the N-body problem has Newton mechanics in the center- 
of-mass frame as non- relativistic limit. At the same time it is the special relativistic limit 
of the rest-frame instant form of canonical metric and tetrad gravity [7], when the Newton 
constant G is turned off. Both in relativistic and non-relativistic theories there is a SO (3) 
left action on the (6iV — 6)-dimensional phase space of the canonical relative variables with 
respect to the 3-center of mass, generated by the non-Abelian Noether constants for the 
angular momentum. Correspondingly, the notion of clustering of the spins of sub-clusters 
do exist at both levels, while the Jacobi clustering of the sub-cluster 3-centers of mass is not 
extendible to special relativity. Previously introduced concepts, like dynamical body frames 
(replacing the standard notion of body frames valid for rigid systems), spin frames, and 
canonical spin bases, proper of the Galilean group-theoretical and Hamiltonian description 
of N-body or deformable systems, are likewise directly extended to special relativity. The 
canonical spin frames and a finite number of dynamical body frames (SO (3) right actions on 
the (6N — 6)-dimensional phase space) can be introduced for every N with well defined non- 
point canonical transformations (when the total angular momentum does not vanish). At 
the non-relativistic level, our treatment generalizes the orientation-shape bundle approach 
ofRef.[l]. 

The main results obtained can be summarized as follows: Our procedure leads to: i) 
The relativistic separation of the center of mass of the isolated system and the characteri- 
zation of all the relevant notions of external and internal 4- and 3- centers of mass, ii) The 
construction of 6N — 6 canonical relative variables with respect to the internal canonical 
3-center of mass for the N-body systems (unlike the non-relativistic case, they are defined 
by a canonical transformation which is point in the momenta in absence of interactions and 
becomes interaction-dependent when interactions are turned on), iii) The simplest construc- 
tion of Dixon multi-poles of an isolated system with respect to a center of motion, naturally 
identified with the internal canonical 3-center of mass. This provides, in particular, the only 
way for introducing relativistic tensors of inertia. 

Such results are intermediate steps in view of future developments. Actually, they furnish: 

iv) The natural theoretical background for a future relativistic theory of orbits for a N-body 
system, by taking into account the fact that in every instant form of relativistic dynamics, 
action-at-a-distance potentials appear both in the Hamiltonian and in the Lorentz boosts. 

v) A new definition of canonical relative variables with respect to the internal canonical 
3-center of mass for every field configuration admitting a collective 4- vector conjugate to 
the conserved field 4-momentum. vi) The analysis of the relevant notions of centers of 
motion and associated Dixon multi-poles for an open sub-system of an isolated system, in a 
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way which can be easily extended to general relativity, vii) The theoretical background for 
a post-Minkowskian approximation of binary systems in general relativity The relativistic 
theory of orbits should provide the relativistic counterparts of the post-Keplerian parameters 
used in the post-Newtonian approximation [56]. viii) The theoretical background for the 
characterization of a relativistic rest-frame micro- canonical mean-field thermodynamics of 
N-body systems with long range interactions (Coulomb or Darwin potentials [6]) as it has 
already been done in Ref.[57] for non-relativistic self-gravitating and rotating systems, ix) 
A theoretical background that could be extended to the weak-field general relativistic N- 
body problem after a suitable regularization of self-energies. Actually, this framework has 
already been extended to charged Klein-Gordon fields interacting with the electro-magnetic 
field [26], Dirac fields [58] and relativistic perfect fluids [27]. Finally: x) Parametrized 
Minkowski theories in non-inertial frames are prepared for a systematic study [9] of the 
allowed conventions for clock synchronizations, the influence of relativistic inertial forces 
and the time-delays (to order 1/c 3 ) for the one-way propagation of light rays (for instance 
between an Earth station and a satellite). 

The extension of our new rotational kinematics to continuous systems like the Klein- 
Gordon field, an extension which should be instrumental to atomic and molecular physics at 
least at the classical level, is under investigation [52]. When applied to the electro- magnetic 
field, these methods could lead to interesting results for the problem of phases [53] in optics 
and laser physics. Also, relativistic perfect fluids have been studied in the rest-frame instant 
form [50] and a future application of the new rotational kinematics might give new insights 
for their description. 

Let us conclude by noting the the non-point nature of the canonical transformations will 
make the quantization more difficult than in the orientation-shape bundle approach, where 
a separation of rotations from vibrations in the Schrodinger equation is reviewed in Ref . [1] . 
The quantizations of the original canonical relative variables and of the canonical spin bases 
will give equivalent quantum theories only if the non-point canonical transformations could 
be unitarily implementable. Up to now, these problems are completely unexplored. 
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